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Preface 



Although Charles Babbage laid down several ideas for computing "en- 
gines," the forerunners of today's computers were largely developed in the 
early 40's as part of the war effort. The Robinson series cryptoanalytic ma- 
chines developed in England in 1941 spawned many families of computers 
still in use today. The MIT differential analyzer and real-time aircraft simula- 
tion project led to the Whirlwind, and eventually to the immensely successful 
DEC family of PDP computers (Programmed Data Processors). And, of 
course, Eniac built at the University of Pennsylvania, was the guiding light 
behind Univac, IBM and many other successful manufacturers. 

Many people today poke fun at these early machines and regard them as 
dinosaur-like relics. However, it is interesting to consider that a large-scale 
computer of about 25 to 30 years ago had about the same amount of power as 
a typical personal computer of today. It was generally not as reliable or user- 
friendly as a personal computer, and, of course, cost tens of thousands times 
as much. Why bring this up? 

Because a computer of the 50's required that the programmer be very 
clever and resourceful to solve problems within the capabilities of the com- 
puter. He did not have vast gobs of memory available, blinding quick calcula- 
tion speed, or random disk access. In other words, he had about the same 
problem to face as you do with your personal computer. 

I do not mean to imply that your personal computer is not a full-fledged 
computer. It certainly is just as much a computer as a room-filling giant of to- 
day. However, because of the relatively small memory, it cannot store a large 
data base. Nor is it suitable for extensive word processing or massive calcula- 
tions. Highly detailed graphics are best left to other machines as well. 

What can we learn from the computing pioneers of the 50's that will help 
us today? Perhaps most important is the discipline of thoroughly analyzing a 
problem, breaking it down into manageable steps, and solving it a step at a 
time. It is also important to determine what can be done "off line" and what 
must be done on the computer. 
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That is what this book is aU about. While it has more than 50 ready-to-run 
programs, the main thing you should look for from the book is an approach 
to solving problems-big and small. Some of the problems demonstrate the 
capability of the computer; others identify its shortcomings. It is important to 
be familiar with both the strengths and weaknesses of your tools so you can 
recognize the types of jobs for which they are suitable (and not suitable). 

This book focuses primarily on mathematical and educational applications 
for the computer. There are many other excellent sources of information 
about other appUcations and for making best use of your personal computer. 
Using the approaches described in this book should enable you to easily con- 
vert programs and use appUcations from other books and magazines such as 
Creative Computing. 

This book is designed to be read with a working computer at hand. While 
there is textual material to be read, the most important things are the experi- 
ments and problems to be tried with your own computer. The book raises 
many questions for which you should try to find answers. There are no an- 
swers to these questions and problems in the back of the book; you should be 
able to discover the answers as you work the problems out on your computer. 

You wUl be able to incorporate many of the routines and approaches in the 
book into programs of your own as you use your computer to deal with "real 
world" problems. Other programs will simply point you in the nght direc- 
tion. And some of the programs in the book are just plain fun. Learn. Expen- 
ment. Have fun! 

Morristown, New Jersey . , tt * t.i 

A * iQoa David H. Am 

August 1983 
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Drill and Practice 



Throughout Hfe, there are certain things that simply must be memorized. 
Obvious things that fall into this category are the addition and multiplication 
tables, the spelling and meaning of words, how to tell time, and the monetary 
system. 

However, depending upon one's chosen profession, there are many other 
things to be memorized. A doctor must know what diseases match what 
symptoms. A chemist must know the gas laws, the properties of elements and 
so on. A pilot must instantaneously know the meaning of readings on scores 
of instruments. 

To memorize a set of facts, you must go over them again and again and 
keep trying different variations. Here is where the computer comes in. It is 
able to present randomly scores of different problems to you for as long as 
you wish. Some programs will automatically adjust to your level of com- 
petence and will grade you; other programs simply present the problems and 
leave the grading up to you. 

There are four programs in this chapter and two in the Science chapter 
which present material in a drill and practice format. Examine the methods 
used in these programs and then make up some drill programs of your own 
for subjects with which you are having trouble, or make up programs for 
other members of your family. 
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Addition Practice 



This program demonstrates a simplified addition drill and practice routine. 
This type of drill is sometimes called computer assisted instruction (CAI), al- 
though CAI can also apply to tutorial and other approaches as well. 

When the program is run, it will first ask "No. of digits?" You enter a 
number and each addend will contain that many or fewer digits. 

The program will present any number of practice problems that you spec- 
ify. The program presents each problem in turn. The program will not pro- 
ceed to the next problem until the current one has been answered correctly. 
After the last problem is answered correctly, the score is printed with an 
appropriate comment. 

There are many improvements and extensions possible in a program like 
this one. For example, you might want to modify the program so it tells the 
user the correct answer after a problem has been answered incorrectly two (or 
three) times. 

A more complicated modification would be to change the program to 
present different kinds of arithmetic problems such as subtraction, multiplica- 
tion, and division. 

Some of these modifications have been made in the next program. 



4 



5 Print chr*<14):rem switch lower cast 

10 Print chr*(147);"Hddition Practi'-*" 

20 inPut "No. of digits" ;n 

30 inPut "No. of Problems" ;p 

40 w=0 

50 x=x+l 

60 q=0 

70 a-int(10'rn*rnda»+I 
80 b-inta0tn«rnd< 1)5+1 
90 Print 

100 Print t«b(8-intao9<aVlo9a0)+l)>;a 
110 Print tab(7-int<lo9<b)/lo9<10>+l));"+";b 

120 Print " " 

130 c»a+b 

140 Print tab<7-int<to9<c)/to9<10)+l)) ; 

130 inPute 

160 if 9»c then 230 

170 q=q+l 

180 if q>l then 200 
190 w=w+l 

200 Print:print "What? Try a9ain, " 
220 9oto 90 

230 Print: Print "Ri9ht! "i 
240 if x>=P then 280 
250 Print "Here's another" 
270 goto 50 

280 Print:Print:Prlnt "Vou got" ;P-w; "correct the first time." 

290 if w<.22«P then 320 

300 Print "Rnd you. missed" J w; 

310 end 

320 Print:print "Oood wor.l"; :end ^TlrViit^f 

No. of Problems? 3 

189 
+ 973 

? 1162 

Right! Here's another 

777 
+ 419 



? 1365 

What? Try again. 

777 
+ 418 



? 1195 

Right! Here's another 

• 030 
+ 810 



? 1S40 
Right! 



Vou 9ot 2 correct the first time, 
find vo'x missed 1 
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Addition Practice, 
Adjusted by Grade Level 

One of the major disadvantages with many drill and practice exercises is 
that they tend to be either boring or frustrating, depending upon the ability of 
the user relative to the level of the material. To compensate for this, a method 
is needed which will adjust the difficulty of the problems to the ability of the 
user. 

Ideally, such a system would weigh the most recent performance most 
heavily but would not ignore previous performance. It should allow a user to 
advance to more difficult problems than his current mastery level. It should 
also continue to give some practice on problems already mastered. 

Some commercial software packages approach these goals along traditional 
lines, i.e., determine in which type of problems a student should receive prac- 
tice by using a compUcated computer managed instruction score recording 
and adjustment system. 

The approach here is more innovative; it uses a single measure for each 
type of problem— call it "estimated grade level"— which meets all of the 
objectives stated above. 

How does it work? The most recent problem presented counts 10% of the 
overall score if it was answered correctly and was over the current user grade 
level, or if it was answered incorrectly and was under the current user grade 
level. Otherwise it is ignored. This may be easier to visualize in the form of a 
chart: 

Answer 
Right Wrong 

Higher than 
grade level 

Problem 

Lower than 
grade level 

At first glance this might look complex and somewhat goofy, however, 
what it really means is that a student is rewarded for doing a problem beyond 
his grade level but he is not penalized if he cannot do it. On the other hand he 
is penalized if he cannot do a problem lower than his grade level, but is not 
rewarded for doing one lower. 

Each problem affects the estimated grade level a little bit, with the most 
recent problems being weighed the most heavily. If the current grade level of 
a student is L and the level of the most recent problem to be averaged in is P, 



Raise student Ignore 
grade level 



Ignore Lower student 



grade level 
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then the averaging formula is simply: 

L = .9L + .IP 

The remaining task before a program can be written is to assign a grade 
level to each problem presented. Unfortunately, this will vary depending 
upon the local school system, the textbook used, and the teaching method. 
Also a huge data base can not be stored in a small computer, so it is desirable 
to devise a simple method of determining grade level for different problems. 
One straightforward approach is to present problems up to one-half a grade 
level over and under where the student currently is. Thus the overall range of 
problems for a student at grade level 3.2 would be 2.7 to 3.7. 

How do we generate the right problems? Consider one type of skill, vertical 
addition. It is normally introduced in the first grade and continues through 
Grade 4 (actually 4.9). The simplest problem in this program is 1 + 1 and the 
most difficult is 999 + 999. Since learning is not a Unear process (it is slow at 
first, and then progresses rapidly), an exponential formula can be used. For 
example: 

Addend = 1.73 x (Grade level)" 

or 



Grade level = V Addend/1.73 

This means that students at various grade levels will be working with the 
following maximum addends: 

Grade level Addend 
1.0 1 
2.0 27 
3.0 140 
4.0 442 
4.9 997 
Now it is a relatively straightforward, although somewhat tedious, matter 
to tie all these elements together in a computer program. 

A few notes about the program. The variable G2 is the problem grade level 
that is always within one-half of a grade level of the current student level, Gl. 
The complicated mess in Statement 290 produces a moving grade level 
average. 

The recording of the grade level and carrying it over to the next lesson is a 
manual process. On a computer system with a permanent mass storage de- 
vice, this would be kept on the system. 

There are many possible changes and extensions to this program. For 
example it could present different types of problems such as horizontal addi- 
tion, vertical and horizontal subtraction, as well as multiplication, division 
and fraction problems. 
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9 Print chr*<14>:reiii »wltch lower case 

20 print chr*<147>;"Hl. To stoP enter 9999 as answer." 

30 Input "Grade lewel";9l 

40 92-«l-.5#rnd<0) 

50 r«lnt<2+1.73#92t4)+l 

60 a*lnta00«92«rnd< 1)5+1 

70 if a>r then 60 

90 Prlnt:print tab(8-lnt<lo9Cb+l)/'lo9<10>+l»;b 
100 Print tab<7-lnt<lo9<a+l)^lo9a0)+l)); + .« 

110 Print " " 

120 r"a+b . 

130 Print tab<7-int(lo9(r)/'lo9(10>+l))' 

140 inPutg 

150 if 9-9999 then 320 
160 if 9-r then 260 
170 u"w+l 

180 if u>l then 210 

190 Print :Print "Uron9, tr» a9ain" 

210 ?rint:?rint "Vou missed that 2 tines. The answer ls";r 
230 w"0 

240 if 92>9l then 300 
250 9oto 290 
260 w-0 

270 Prlnt:Prlnt "Correct!"; 

280 if92<91 then 300 

290 9l-.9#91+. l#s<lr<sqri;r/'2/'1.73)) 

300 Prlnt^Prlnt "Here's another. .." :Prlnt 

310 9oto 40 

320 Print chr«<147) 

330 Print "Okan. So lon9." 

340 9l«<int<100*91>)/'100 

350 Print "New 9rade level";9l 



Hi. To stop enter 9999 as answer. 
Grade level? 4 

349 
+ 73 



7 422 
Correct I 

280 
+ 125 

? 405 

Correct ! 



296 
+ 55 

? 341 
IJron9< trK aOaln 



296 
+ 55 



? 351 
Correct I 



200 
+ 177 

? 377 

Correct ! 



76 
-f 211 

? 9999 

Okav. So lon9. 

New 9rade level 3.93 
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"You won't believe this, but Zarg can actmlly add without a computer. 
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Time/Speed/Distance Problems 



As well as being able to present simple numerical drill and practice prob- 
lems, the computer can present word problems for solution as well. 

In this program the formula relating time, speed and distance was applied 
to a problem involving both a car and train. 

The problem can be stated as follows. A car traveling C miles per hour 
(computer generates an integer 40 through 65) can make a certain trip in D 
hours (computer generates an integer 5 through 20) less than a train traveling 
at T mph (computer generates an integer 20 through 39). How long does the 
trip take by car? When the two simultaneous equations are solved they pro- 
duce the single equation for the answer shown in Line 60. 

Notice the calculation in Line 70. This calculates the percent difference be- 
tween the actual answer and the one entered by the user. 

Notice also that the computer calculates the correct answer in Line 60 and 
prints it (on the screen) in Line 130. This answer may have many decimal 
places; as the program is written it is rounded off to two decimal places. If the 
0.5 was not added in Line 60, the number would be truncated and not 
rounded off. This is an important calculation and one which you will find in 
many other programs throughout the book. 

Consider other problems that can be used as the basis for this kind of drill 
and practice exercise. Teachers, for example, might wish to have students 
write drill and practice programs on their own. Different problems could be 
given to one or a small group of students to serve as the basis for a program. 

After the programs are written, students can try out the programs of other 
class members. This approach ensures that each student not only understands 
the type of problem assigned to him, but also gets practice in solving other 
problem types as well. This is an effective technique for stimulating interest as 
well as for learning how to solve word problems. 
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5 Print chr*<14):rtm switch tower cut 
10 c-lflt<rnd(l)#25)+40 
20 d-lnt<rnd<l)*15)+S 
30 t-lnt<rnd<l>*19)+20 
35 Print chr*<147>; 

40 Print "fl c*r »t";cJ"iiiPh m«ke» a triP ln";d;"hrs iMS "; 

45 Print "th«n * trtln jit";t;"nPh 

50 Print: input "C«r trlP hr»";* 

60 v-int(.5+100*<d*t/'Cc-t)>>^100 

70 e'«lnt<«b»<<v-«)«100^«)+.5) 

80 Print 

90 If e>5 then 120 

100 Print "Good. Vou were within". "Percent. " 
110 9oto 130 

120 print "Sorrs. Vou were off bs";*; "Percent. " 

130 Print "Correct time"" J v 

140 Print^inPut "Another <a,n)";rl 

150 If z*="s" or 2*«"V" then 10 

160 Print:print "Oku. So lon9." 

170 end 



fl c*r at 53 mPh makes a trip in 5 hrs less 
than a train at 23 mPh 
Car trip hrs? 4 

Good. Vou were within 4 Percent. 
Correct tlm»« 3,83 

Another <i(,n)7 v 

fl car at 58 nPh makes a trlP in 14 hr* less 
than a train at 31 mPh 
Car trip hrs? 13 

Sorry. Vou were off bv 24 Percent. 
Correct time" 16.07 

Another («,n>7 y 

fl car at 56 mPh makes a triP in 11 hrs less 
than a train at 26 mPh 
Car trip hrs? 7.25 

Sorry. Vou were off by 31 Percent. 
Correct time" 9.53 

flnother <yjn)? n 

Okay. So lonS. 
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Kinematics Problems 



Kinematics relates to the dynamics of motion of bodies apart from consid- 
erations of mass and force. Like many other types of problems, these can be 
generated and presented by the computer to provide practice m solvmg them. 

This program presents a simple kinematics problem for solution. The com- 
puter generates a new value for each problem. The problem is as follows. A 
ball (or any other object) is thrown up at velocity V meters per second (com- 
puter generates an integer between 5 and 40). The user must then calculate 
three factors about the resulting flight of the ball: maximum height, time until 
it returns, and velocity after T seconds (computer generates a time less than 

the total flight time). . .u- . 

The key benefit in using a computer to present problems of this type is 
motivation. The calculations required of the user are no different than those 
in the back of a chapter in a book or those on homework assignments. How- 
ever, answering them when they are presented by the computer seems to 
make it more of a challenge and, frankly, more fun. 

The computer program checks each of your responses to see if it is within 
15% of the correct answer. If it is, your answer is considered correct. You 
may wish to change this percentage to require more or less accurate calcula- 
tions. You may also wish to change the computer calculations to round off to 
one or two decimal places. 




5 Print chr*<14>:rtB twitch loutr case 
10 0-e 

20 v»6+lnt<35«rtid<l)) 

30 Print chr*<147);"B«U i» thrown uP jit" ; v; "mettr»/'»ec" 
40 »".05«v12 

50 input "Wh«t 1* tht maxlnuM hti9ht";9 
60 9o«ub 210 
70 ««v/5 

80 Print: inPut "Whtti will it come down <»ec)";9 
90 gosub 210 

100 t"2+lnt<2#v#rnd<l>>/'10 
110 if t>- v/S thtn 100 
120 fv-iemt 

130 Print :print "Wh«t is velocltK *t";t;"secs"; 
140 inPut9 
130 9osub 210 

160 Print:print<i;"ri9ht out of 3" 
170 if q>l then Print "Not b*d." 
180 Print :inPut "flnothtr <K,n)";z* 
190 if z«-"K" or z«»"V" then 10 
200 end 

210 if «b»<<9-a>/'a) <.15 then 240 

220 Print "Not even close..." 

230 9oto 260 

240 Print "Pretts* 9ood. " 

250 q«q+l 

260 Print "Correct answer •";« 
270 return 



B«U is thrown uP «t 24 meters/sec 
Uhat is the max i nun hei9ht? 30 
Pretty 9ood. 
Correct answer •28.8 

When will it cone down (sec)? S 
Pretts* 9ood. 
Correct answer "4.8 

What is ueloclts* at 2.3 sees? 2.3 
Not even close. . . 
Correct answer » 1 

2 riSht out of 3 
Not bad. 



Ball is thrown uP at 13 neters/sec 
What is the max i nun hel9ht? 8 
Pretty 9ood. 
Corr-ect answer « 8.45 

When will it come down <sec)? 2.5 
Pretty 9ood. 
Correct answer "2.6 

What is velocity at 2.3 secs?-9 
Pretty 9ood. 
Correct answer "-10 

3 ri9ht out of 3 
Not bad. 



Another (y^n)? y Another i;y<n)? n 
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2 

Problem Solving 



In many courses in school, the textbooks present a great variety of devices 
for solving the problems that have been neatly grouped together at the end of 
each chapter. Typically these devices consist of formulae, equations, rules and 
theorems. After a careful study of these devices, teachers give exams which 
test your ability to recall them. 

But what do you do if you are faced with the more realistic situation of not 
being told what device is likely to solve which problem or, worse yet, of hav- 
ing forgotten how to use a technique altogether? Is all hope lost? Of course 
not, although some people seem to believe that it is. 

In this chapter, several of these nasty devices mentioned above are pre- 
sented. For example, there are devices written into computer programs that 
will solve a quadratic or exponential equation and others that will calculate 
the roots and draw a plot of any function. 

However, you must remember that while these devices are useful in solving 
certain problems, the really important thing is to understand the underlying 
logic and approach. Then when it comes time to solve real world problems 
you will be better prepared to face them. Incidentally, many of the methods 
and devices presented in this chapter are used in later chapters to solve other 
kinds of problems. 
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How Many Tickets? 



Here is a problem. At a school raffle to raise money, the organizers have as 
a prize an electronic game for which they paid $18.00. To add interest to the 
raffle, the organizers have decided to sell tickets for an amount (m cents) 
equal to the number on the ticket for tickets numbered 1 to 50. For ticket 
numbers over 50, the price is 25 cents each. The organizers want to know 
how many tickets they must sell to exactly break even. 

It is probably easiest to visualize a problem of this sort with a flowchart. In 
the flowchart, T will equal the total money collected and will increase as 
more tickets are sold. The ticket number is N. When T equals or exceeds 
$18.00, N will be the answer. 

Note that the flowchart has two logical branching points (if statements in 
the program). The first compares the current ticket number to 50; if it is less, 
the ticket number is added to the total whereas if it is greater than 50, the 
total is increased by 25 cents. 



T=0 




The second branch point compares the total amount collected, T, to $18.00 
(actually 1800 cents). If T is equal to or greater than 1800, the break even 
point has been reached and the values of N and T are printed (on the screen). 

A problem of this kind can be done by hand, however, because of the repet- 
itive additions it is quite tedious. Also, doing it by hand frequently leads to an 
answer of 72 rather than the correct answer of 71. Try it yourself and see 
what you get. 

10 t-0 

20 for n"l to 2000 
30 if Ti<«50 then 60 
40 t-t+25 
30 9oto 70 
60 t-t+n 

70 if t>1800 then 90 
80 next n 

90 Print n; "TICKETS SOLD" 
100 Print "COLLECTED «";t^IB0 



71 TICKETS SOLD 
COLLECTED f 18 



Many problems can be solved quickly and correctly with a computer using 
logical analysis and a flowchart. More complex problems may have to be bro- 
ken down into additional steps and require a longer flowchart, but the ap- 
proach is fundamentally the same. 

Here are two problems for you to solve. 

The diameter of a long-playing record is 12 inches. The unused center has a 
diameter of 4 inches and there is a smooth outer edge 1/2 inch wide around 
the recording. If there are 91 grooves to the inch, how far does the needle 
move during the actual playing of the recording? 

A movie theater charges $2.50 for an adult admission and $1.00 for a child. 
At closing, the cashier counted 385 ticket stubs and had $626.50 in cash. 
How many children entered? 
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Drinking and High Blood Pressure 



This program illustrates how several simple equations can be put in a com- 
puter program to solve a more difficult overall problem. 

Here is the problem. In a survey of 1000 adults, it was found that 35 had 
high blood pressure. Of those with high blood pressure, 80% drink 15 oz. or 
more of alcohol per week. Of those without high blood pressure, 60% drink a 
similar amount. What percent of drinkers and non-drinkers have high blood 
pressure? 

This problem requires the solution of several simple equations. They could 
all be combined into one large equation, but it may be easier to understand 
the approach (and change variables later on) by writing a program with five 
separate equations. 

If H equals the number of people with high blood pressure, then H = 35 
(Line 10). Letting HI equal the number of people with high blood pressure 
who drink leads to HI = .8 x H (Line 20). 

Letting LI equal the number of people with low blood pressure who drmk 
yields LI = .6 x (1000 - H). Then, the total number of drinkers, D = HI + 
LI. 

Finally, the percentage of drinkers with high blood pressure is X = HI x 
100 / D. 

The program solves the problem in a jiffy. The solution for this type of 
problem can be easily written directly in Basic without any need for a flow- 
chart or detailed analysis. Recognizing this type of problem readily will save 
a great deal of pencil pushing time. 



10 h»35 

20 hl«.8*h 

30 U«.6*i:i080-h) 

40 cl«hl+U 

30 x-hftlOe/d 

€0 xi'h«100/'i000 

70 Print "HIOH BLOOD PRESSURE" 

80 Print "POPULATION" ;xl;"Ji" 

90 Print "BRINKERS ";x;"!<" 



HIGH BLOOD PRESSURE 
POPULATION 3.3 5i 
DRINKERS 3.76606261 X 
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Here is a problem that doesn't require a single equation but makes use of 
the approaches discussed so far in this chapter. Can you solve it with three 
Basic statements? 

In early January, a shopkeeper marked down some calendars from $2.00 to 
a lower price. He sold his entire stock in one day for $603.77. How many did 
he have? 

Here's another easy one. A town in India has a population of 20,000 peo- 
ple. Five percent of them are one-legged and half of the others go barefoot. 
How many sandals are worn in the town? 




Two Simultaneous Equations 

So far in this chapter, only problems with linear equations have been 
considered. But the computer can be used to solve much more difficult equa- 
tions. In fact, it is in problems involving second and third degree equations, 
exponentials, and the like where the computer really starts to pay off. Con- 
sider the following two simultaneous equations: 

2X = 16y 3X = 27y 

It is not at all easy to solve these two equations by hand. But a simple Basic 
program can be written to solve the equations using trial and error. This is 
sometimes referred to as a brute force approach because every possible 
combination of numbers between an upper and lower limit is tried until a 
solution is reached or until the program runs out of values. 

5 Print "SOLVES TWO SIMULTRNEOUS EQUATIONS" 
10 for x"l to 4 
13 for x«l to 4 

30 if'"21x<>16*n^3 or lflt<3tx)<>27«» thtn 60 
40 print "x"";x/'K-";s< 
50 »toP 
£0 ntxt n 
70 ntxt X 

80 Print "NO INTEGER SOLUTION" 



SOLVES TWO SIMULTANEOUS EQUATIONS 
X* 4 «" 3 



In this particular case, a solution is reached rather quickly with x = 4 and 
y = 3. However, if in the second equation the y coefficient is changed slightly 
from 27 to 28, the computer will try 10,000 possible solutions before finally 
concluding that no integer solution exists— at least within the range of to 
100. Warning: this will run for a very long time. 

30 if2tx01S#»/3 or lnt<3tx)028«i) then 60 



SOLVES TWO SIMULTANEOUS EQUATIONS 
NO INTEGER SOLUTION 



20 



Although the trial and error (brute force) approach is widely used, it is 
highly inefficient. In general, a systematic or guided trial and error approach 
is preferable to one that simply tries every possible solution. However, for 
some problems the simple "try every value" approach may be appropriate. (A 
comprehensive discussion of trial and error approaches can be found on pp 
36-40 of Computers in Mathematics: A Sourcebook of Ideas.) 

Three problem solving approaches have been discussed so far. Remember- 
ing them, how would you do this problem? A boy and his sister visited a farm 
where they saw a pen filled with pigs and chickens. When they returned 
home, the boy observed that there were 18 animals in all, and his sister re- 
ported that she had counted a total of 50 legs. How many pigs were there in 
the pen? 




Quadratic Equation Solver 



For any values A, B, and C of a first degree quadratic equation (Ax^ + Bx 
+ C = 0), this program will compute the roots of the equation. The solution 
is based on the quadratic theorem which solves for roots with the following 
formula: 

_ -B+VBMAC 
^ 2A 

Assuming that A, B, and C are real numbers, the following principles 
apply: 

1. If B^ - 4AC is positive, then the roots are real and unequal. 

2. If B^ - 4AC equals 0, then the roots are real and equal. 

3. If B^ - 4AC is negative, then the roots are imaginary and unequal. 
The program takes into account all these possibilities and correctly identi- 
fies the type of roots along with their values for any set of coefficients. 

Is this program useful by itself? Except for solving quadratic equations for 
algebra class, probably not. However, as a routine in a larger program to 
solve quadratic equations that might be encountered, it could be very useful. 



10 (>rlnt: Print "QUflDRflTIC SOLVER" 

20 InPut "fljB,C";»<b/C 

30 r«bt2-4*««c 

40 If «O0 th«n S0 

30 Print "FIRST DEGREE EQUATION" : «oto 150 
ee Print "ROOTS ARE "; 
70 if r<0 thtn 120 
80 if r-0 then 118 

90 Print C-b+»<lr<r))/'2««; (-b-*<Ir<r> V2«» 
100 Soto 150 

110 Print -b/'2*i!floto 150 

120 Print "IMflOINHRV" 

130 Print <-b/2*«);"+";s<lr<-r)/2«J:;"*i" 

140 Print <-b^2#t>; "-";sflr<-r),'2ilHi; "#i" 

130 input "ANOTHER <V<N>";»« 

160 if *«-"«" thtn 10 

170 »toP 



QUABRflTIC SOLVER 
A>B.C7 1,2.1 
ROOTS ARE -1 
ANOTHER <V.N>? V 

QUADRATIC SOLVER 
A,B,C7 0,2,4 
FIRST DEOREE EQUATION 
ANOTHER <V,N)? V 



QUADRATIC SOLVER 
fl,B-C? 2,5,2 
ROOTS ARE -2 -8 
ANOTHER (V,N)? V 

QUADRATIC SOLVER 
A,B,C? 4 , 2 , 4 
ROOTS ARE IMAQINARV 
-4 + 15.4919334 #1 
-4 - 15.4919334 #1 
ANOTHER <V,N>? N 
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Exponential Equation Solver 

Another general routine for solving a particular type of equation is this one 
to solve for an exponent in an exponential equation. 

Given the values A. B, m. and n. this program will solve for x in any 
exponential equation of the form: 

For example, the program will solve any of the following problems: 

1. = 40 

2. = 7.6 

3. 17X-3 = 8.12 

4. 11 '-^x = 247 

If you were to solve an exponential equation by hand, you would probably 
go through the following steps: 
5^ = 40 
log -- log 40 
X log 5 = log 40 
X = log 40/log 5 
X = 1.6021/.6990 = 2.292 
However, in more generalized form, the solution for x is: 

log B 
X = 'og A 



M-(N/M) 



The data for the four problems listed above were entered into the four data 
Statements (100-130). 

The program as it is presented here can be improved in several ways. First, 
it always solves the same four equations. How can you generalize it to solve 
for other equations? Second, if you were to make use of this routine in an- 
other program, you would probably not be able to use a read statement; how 
could you get rid of it?. 
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5 Print chr«(14>:i"ti«i »ultch lower c»m 

le Print chr«<147>;"exP0tnnti«l Eflu*tlof> Solvtr" 

20 Print 

30 Print " fl B M N X" 

40 Print " " 

50 rt«d »,b/n/n 

60 x"<lo9(bVto9(»)>/'iii-<in<'i»> 

70 x«lnt<.5+l00«x)/'100 

80 Print «;t«b<7);b;t*b<14>;i»i;t*b<21);n;t*b<27);x 

90 Ooto50 

100 d«t« S<40,1.0 

110 d«t* 5-7.6<3/l 

120 d«t« 17,8.12.1.-3 

130 d«t« 11,247.-2.1 



ExPon»ntiil E<lu«tlon Solver 



fl 


B 


M 


N 


X 


s 


40 


1 





2.29 


5 


7.6 


3 


1 


.09 


17 


8.12 


1 


-3 


3.74 


11 


247 


-2 


1 


-.65 
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Roots of Any Function 



This program will find the roots of a function, any function! The function 
may be linear, quadratic, cubic, trigonometric or any combination as long at 
it can be represented in the Basic language. The program as it appears here 
finds the roots between -20 and 20 although you can change these boundaries 
in Statement 50. 

The method used involves evaluating the function at small incremental 
intervals, finding places where the value of the function changes sign and 
then, by successive approximations, finding the zero point. This approach 
borrows from Newton's method in the final narrowing down but, unlike 
Newton's method, will not fail to converge in the event one makes an unlucky 
first guess. 

Before running the program you must first type in your function in state- 
ment 30. For example, 

DEF FNA (X) = 2*XT3+11*XT2-31*X-180 

DEF FNA (X) = X-4 

DEF FNA (X) =SIN(X)-.5 

You may have to refer to the Basic manual with your system to see exactly 
how a function should be stated. 

The routine used in this program is very powerful and could possibly be 
used as a subroutine in many other programs. 

How can you use this program to help you solve this problem for x? 
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5 Print chr*a4):rem twitch lowti" c»«e 

10 Print chr*( 147); "Root* of imi function. Define function in line 30 

30 dtf fn«<x)«2*xt3+U«xt2-31*x-180 

40 zl>-200 

30 for 1«-19.9 to 20 

£0 if »9n<fn*(l>)»«9n<fnt{i+n) thtn 230 

70 k«l 

80 J-i'^1 

90 if fn*<kKfn»<J) then 130 

100 z-k 

110 k-J 

120 J»2 

130 z-<k+J5^2 

140 if fnji<z><0 then 170 

1S0 J-z 

160 Soto 180 

170 k-z 

180 if «bsCfn«<z>». 00005 then 130 

190 2"s9n<z)«int<*b*<z>*10000+. 85>/'10000 

200 If z>zl then 230 

210 Print "f(";z;")«0" 

220 zl-z 

230 next i 



Roots of *ns function. Define function in line 30 
f<-5 >-0 
f<-4.5 -0 
f< 4 )-0 



30 def fn«<x)"x-4 

Root* of tnK function. Define function In line 30 
f< 4 )-0 



def fn*<x)»»in<x)-,5 



Roots of «n» function. Define function in line 30 

f<-18.3259 >-0 

f<-16.2313 >-0 

f (-12.0428 )-0 

f<-9.9483 )-0 

f (-5. 7596 )»0 

f<-3.6651 >«0 

f< .5235 )»0 

f( 2.6179 )»0 

f< 6.8067 )"0 

f( 8.9011 )'»0 

f< 13.0899 5-0 

f( 15.1843 )-0 

f< 19.3731 )"0 



27 



Plot Any Function 



Here is a nifty program that will produce a plot of any function on the 

Te^rl rTntg th.s program, you must type in your func.on in lin^m 
T ike the orevious program, which finds the roots of any function, this one 
tmV'ot aTyTnctL^You' must tell the program between w"^^^^^^^^ 
want the function plotted, i.e.. a minimum and maximum value of the x co 
ordinate. You also input the x plotting increment you wish 

AS it appears here, the program does not allow the user to select the y co- 
ordinates; the program plots y values between -30 and 30. 

Here are several functions you might want to try plotting: 



X Limits 
-15 15 



■5 
-5 



5 
6 



-30 30 
-5 6 



Increment 
1 

. 25 
. 5 
1 . 5 
1 



Function 

DEF FNA(X) = 2«X 
DEF FNA(X) = 30*SIN(X) 
DEF FNA(X) = X-X*X 
DEF FNA(X) = 30«EXP(-X*X/100) 
DEF FNA(X) = X*X-X 
DEF FNA(X) = XT2-X-15 

The last function listed is the one plotted in the sample run with the pro- 
gram. Exponential functions are a great deal of fun and sometimes lead to 
unexpected and interesting results, particularly when combined with trigo- 
nometric functions. Experiment! Have fun! 



5 Prifltchr*<14):rei« switch lowtr «»» f„„rfioY> it, \ir<t SB." 

10 print chr«( 147); "Plot »n« fumctioti. Define functlom im 

20 print 

30 def fn»<x>«xt2-15 

40 inPut "InltUt v«lue of x";xl 

50 input "Fln»l vjilue of x";x2 

S0 inPut "V»l of »teP";» 

70 Print t«b<19);"X«";xl 

80 for X"xl to x2 steP -s 

90 If *bt<x><. 00001 then 160 

100 v>fni<x)'«'20 

110 if s»<40 then 130 

120 y>39 

130 if K>20 then 200 

140 print t»b<i«>;"«";t*b<20V':" 

}l0 Prfnt'"?V"rfor t-l to 36:Print "-".-..Kt t^Print "*v" 
170 Print " -20 -10 *10 *^ 

180 x"x-» 
190 9oto 210 

200 Print t»b<20);":";ttb<»);"#" 
210 next X 

220 print t»b(l9>/ "X-";x2 
230 9oto230 
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Pl«t jpny function. Define function in 
1 iiie 38 , 

Initial valiM of x? 6 
Final value of x? -6 
Ual of step? 1 

X= 6 



-V — 

-za 



-le 



— +v 

+26 



x=-e 



Pl4>t any function. Define function in 
liAC 39. 

Initial salue of x? 6 
Final value of x? -6 
Ual of step? ^5 



X> 6 



« 

-V 

» 



+v 



X=-S 
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Sets and Repetitive Trials 



For solving relatively simple problems, the computer may not be any help 
at all. In fact, it may take more time to write a program to solve a problem 
than it would to solve it by hand or with a calculator. This chapter should 
help you recognize problems that are suitable for computer solution and 
those that are not. 

In some problems, you may think that the computer will not be any help. 
However, one thing writing a computer program will always do is force you 
to reason out the approach to solving the problem logically and precisely. The 
computer can't solve problems unless it is told exactly how to proceed; hence 
you must understand a problem completely before you can program it for the 
computer. 

Several of the sections in this chapter discuss sets of data. While the sets 
used in the examples have relatively few elements or values, you should bear 
in mind that real world problems often have thousands or millions of pieces 
of data, and the only practical way to solve problems of this size is with a 
computer. For example, consider how you would most efficiently schedule 
the shipments on a railroad train leaving Boston with 4000 diverse cargoes 
bound for Phoenix and 780 points in between. Now consider that there are 
200 freight trains per day leaving Boston. Now add to that the 10,000 other 
trains leaving other cities every day that must use the same network of track 
and you can see that dealing with real world sets of data is no easy task. 
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Group of Girls and Boys 

In the previous chapter we said that there may be some problems for which 
brute force trial and error is appropriate. This would be the case if the prob- 
lems were relatively small and trying every possible solution would not tie up 
a great deal of valuable computer time. Here is a problem involving two lin- 
ear equations that lends itself to a trial and error approach. 

The problem is as follows: When 15 girls leave a group of boys and girls, 
there are two boys for every girl (lucky girls). Next, 45 boys decide to leave; 
then there are 5 girls for every boy (lucky boys!). How many girls were there 
in the group before anyone left? 

Before rushing to the computer, you must recognize that this problem re- 
quires the solution of two simultaneous equations. If G equals the original 
number of girls and B the original number of boys, then the two equations 
are: 

(G - 15) X 2 = B 
(B-45) X 5 = (G- 15) 



5 print chr«<14>;i"tiii iwitch lowtr c*»e 
28 Print" Pro9rtn Is 3L0W< bt Putltnt " 

30 for «■! to lee 

40 for boi to 100 
S0 1"1+1 

60 If 2#<9-15>Ob thtn 100 
70 If 5«<b-45K><9-15> th»n 100 
80 print 9;"01rl» ";b;"BoK»" 
90 eoto 130 
100 ntxt b 

110 next 9 .-,.,»«„ 
120 Print "No lnte9er solutloti from to 100. 
130 print l;"Trl»l»." 



ProSrim i« SLOW- be Pttitnt. . . 
40 OlrU 50 Bous 
39S0 TrUlt. 



The computer program uses two FOR loops (Statements 30-110 and 4( 
100) to try every combination of values for B and G between 1 and 100 until 
solution is found or until the program runs out of values. The variable 
(Statement 50) is a counter which records the number of trials required 1 
reach a solution. 
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The program is straightforward and finds a solution after 3950 trials. How- 
ever, it would have been a simple matter to substitute the value of B from the 
first equation in the second one and quickly solve the problem by hand or 
with the aid of a calculator. It is important to recognize that if a problem can 
easily be solved by other methods, the computer offers little or no advantage. 

Try this problem. You may or may not want to use your computer. If Mat- 
thew can beat Jeff by one-tenth of a mile in a two-mile race and Jeff can beat 
Steven by one-fifth of a mile in a two-mile race, by what distance could Mat- 
thew beat Steven in a two-mile race? (Hint: the answer is not 3/10 mile.) 
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Brown's Books 



The use of a trial and error approach can generally be improved signifi- 
cantly if the combinations to be tried can be narrowed down in some way. 
The solution to this problem illustrates how the speed of obtaining a solution 
can be improved well over 100 fold by combining equations and eliminating 
certain solution possibilities. 

Here is the problem. Brown sold 48 books at a flea market, some for $3 
each, some for $5 each and others for $8. He collected a total of $175. He 
remembered having an even number of $5 books. Can you determine how 
many of each kind of book he had? 

The equations for solution are (letting T equal the number of $3 books, F 
the number of $5 books, and E the number of $8 books): 

T + F + E = 48 

3*T + 5*F + 8*E = 175 

The first program was written simply to try all possible combinations of T, 
F, and E from 1 to 48. It yields three solutions for the problem, although the 
two solutions with an odd number of $5 books can be eliminated leaving just 
the one desired solution. 

5 Print chr*<14>:r»ni iwltch lower c*lt" 

10 Print "Broum'» Books." 

20 Print "L-O-N-0 run tint" 

40 Print "*3 «5 «8" 

30 for t"l to 48 

g0 for f"l to 48 

70 for t«l to 48 

80 If <t+f+e)048 thtn 110 

90 if t*3*f«3«e«80175 thtn 110 

100 Print t;t«b<4>;f;t»b<10);t 

110 ntxt c 

120 n«xt f 

Brown's Book*. 
L-O-N-0 run time 



«3 


*5 


$9 


40 


3 


5 


37 


8 


3 


34 


13 


1 



This program took approximately 22 minutes and 13 seconds to run on the 
Commodore 64 computer. A typical minicomputer (PDP-8/e) could run this 
problem in about 7.3 seconds. In either case, this is a long time to tie up the 
computer. 
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It is rather easy to combine the two equations into one by solving for T. 
The single equation is then: 

2*F +5*E = 31 

In this equation, the limits can be reduced (from 48 used in the first run) 
since F cannot possibly be greater than 31/2 or 15.5 and E cannot be greater 
than 31/5 or 6.2. Making the appropriate program modifications leads to the 
second program. 

20 

se 

50 for f-1 to 15 

70 for t"l to 6 

80 If <f«2+»#5)09I then 110 

90 t«48-f-t 

130 

Broun 't Books. 



L-O-N-0 run 


tin* 


«3 


«5 


*8 


40 


3 


5 


37 


8 


3 


34 


13 


1 



Using this program produces a dramatic improvement in the time to solu- 
tion. On the Commodore 64, the time is approximately 2.5 seconds and on 
the PDP-8, about 0.16 seconds. 

Since the problem states that F must be even, a final modification which 
steps F by two in Statement 20, can be made. This version of the program 
takes only 1.25 seconds to run on the Commodore 64 and 0.06 seconds to run 
on the PDP-8. 



30 Print "F Even" 

60 for f«2 to 14 »ttP 2 

Brown's Book*. 
F Ewtn 

*3 «3 tB 

37 8 3 



Notice the enormous improvement in computing time required for a solu- 
tion, over 1000 fold on the Commodore 64 and 100 fold on a PDP-8. Brute 
force certainly is inefficient! It is generally worthwhile to think through most 
problems, particularly big ones, before rushing to the computer. The com- 
puter may be fast, but we just improved its performance by 1000 times by 
using a little common sense. 
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Intersection of Sets 



Two sets of numbers can be combined to yield a third set by the operation 
of intersection. The intersection of two sets A and B is the set that contains all 
elements that belong to both A and B. It does not contain any other elements. 
The intersection is usually written A n B. . » - „ ^ 

For example if M = < 0,2,4,6 > and K - < 1,2,3,4 > , then M fl K = 

< 2 4 > 

This program finds the intersection of two sets of numbers. It has been 
written to find the intersection of the two repetitive sets described m State- 
ments 30 and 40. In the sample run, Statement 30 describes the set 
X = < 1,3,5, ... 19 > and Statement 40 describes the set 
y = < 2,5,8, . . . 29 > . 
Notice that successive values of x increase by 2 and y by 3. 



10 print chr«<14>:ren switch lower c»»e 

28 Print "S»t» x< ii int«r»ect *t Point* • . 

30 for X"l to 19 »teP 2 

40 for »>2 to 29 sttP 3 

50 if thtn 90 

£0 next V 

70 next X 

80 end 

90 Print x; 

100 Qoto 70 



Set* X, V Intersect tt point*: 5 u 17 



However, if the set cannot be so neatly described, it may de^'^^^le to 
rewrite the program to examine any set of data. This is done with the read 
statement which reads into x the data points in the data fte"»«"t.The pro- 
gram is set up to use the same y set as the first program, but the x set is de- 
fined in the data statement. ^ . * f ic a 

You should be able to see from the first combination of sets hat if there is a 
numerical pattern in the sets which intersect, then there is also a Pattem 
the resulting intersecting set. In the example, the x values mcrease by 2 and 
the y values by 3, hence the values in the intersecting set increase by 2 x 3 - 
6 Although the intersection of these sets could easily be calculated by hand, 
the computer can be an aid in evaluating more complicated sets. 
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5 Print chr«a4):reM switch lowtr c*»e 
10 Print "Set* Xi k intersect at Point*:"; 
20 read x 

30 for «"2 to 29 steP 3 

40 If x*« then 80 

50 next v 

S0 9oto 20 

70 end 

80 Print x; 

90 Soto 60 

100 data 2, 3, 8, 9, 14, 15, 20- 221 < 26/ 27 



Set* X, K Intersect at Points : 2 8 14 20 26 
?out of data error in 20 
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Prime Factors 

A prime factor is a positive integer that has no factor except itself and one^ 
The first ten prime factors (or numbers) are 2, 3, 5. 7, 11,13 17, 1% Z3 ano 
29. The definition gives the basic method for determmmg ^'hether a number 
is prime: divide by all smaller integers down to 2, testmg whether the remam- 
der is zero for at least one of them. If not, the number is pnme 

But this is highly inefficient. It is obvious that a number is not prme if t is 
any even number greater than 2; hence only odd divisors need ^o hcined^ 
Also, it is not necessary to try divisors greater than the square root of the 

" W the division method is inefficient, various schemes have been devised 
to avoid division. The basic idea underlying all such -he-es is call^^^^^^^^^ 
sieve of Eratosthenes (276 B.C.-195 B.C.). Imagine a list of odd numbers 
from 3 up. Strike out every third number after 3, every fifth number after 5, 
and so on. This will leave only prime numbers. 



R9«ln..."i4oto 50 



5 print chr»<14>:i-tm twitch louti- c»it 

le dim ji<iee> 

20 c«0 

30 print "Prlmt Ftctor*" 
40 Print »Ent«r to »toP" 
50 Infut "Vour numbtr'Sn 
50 If n*0 thtn ttoP 
70 n"n 
90 x"0 

90 if n>0 thtn 110 
100 Print "Invalid. 
110 i-1 
120 i-1+1 

130 If l>i« thtn 190 
140 if i«/'i01nt<m/'i) thtn 120 
150 x"x+l 
160 »<x>-i 
170 

100 eoto 140 
190 if x-1 thtn 240 
200 for l"l to X 
210 Print *<t>; 
220 ntxt I 
230 Soto 250 
240 Print n;"i» Print. "J 
260 Print 
270 9oto 50 
280 tnd 



Print Fjictort 
Enttr to itoP 
Vour nunbtr 105 

3 5 7 
Vour nunbtr 72 

2 2 2 3 3 
Vour nunbtr 362 

2 181 

Vour nunbtr 09 
89 U Print. 
Vour nunbtr 
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The program here finds the prime factors of any integer, or prints out ".N is 
prime" if the integer has no proper divisors. 

Run this program for a large number of different integers and see if you 
can discover relationships between numbers and their prime factors. You 
should also try to figure out the method employed in the program to find the 
prime factors of any integer. To do this, you might want to draw a flowchart 
to show what is happening in the program. This will help you see the method 
used to find a prime factor and might help you in writing a program to gen- 
erate primes. 

In writing a program to generate prime factors, you can use the sieve 
method. However, as the numbers become very large, you will have to figure 
out a way to represent integers with more digits than your computer can han- 
dle at one time. (One approach is described on pp. 19-21 of Computers in 
Mathematics.) 

Goldbach was a mathematician who made a conjecture that every even 
number greater than 4 can be written as the sum of two prime numbers (16 = 
11 + 5, 30 = 17 + 13, etc.). No one has ever proved it but no one has dis- 
proved it either. That is why it is called a conjecture. Can you write a pro- 
gram that will prove or disprove this conjecture? Or how about writing a 
program to prove Goldbach's conjecture for even numbers up to 50? You 
should be able to write this program with 12 or fewer statements. 

Here is another problem involving prime numbers. Assume a life span of 
80 years. In what year of the 20th century (1900-1999) would a person have 
to be born to have the maximum number of birthdays occurring in prime 
years? The minimum number? 
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Greatest Common Divisor 

The greatest common divisor of a set of numbers is, as its n^'"; ^P^^^^^^^^^ 
greatesUnteger that will divide into a set of two or -^^^^^^^ d^^ oM 
pie, the set of numbers 12, 20, and 28 have a greatest common dmsor 
Nothing larger than 4 will divide evenly into all three numbers 

This program will find the greatest common divisor for -^y ^l^'^^m in 
To run H, you simply input the number of integers - 
when requested and let the program calculate the GCD. The heart 
calculation is in Statement 150. „f nnmhers'? Can vou 

Do you know the meaning of a relatively pnme ^J'^""^^^^^^^^^ 
figure out the meaning from the third sample run of the P'^^^'^^l''^ 
runs of vour own"? How is a set of relatively prime numbers d fferent from a 
setXre fLtors? Can you find a set of 10 integers that xs relatively pnme? 



5 Print chi"«<14>:reiii switch lowtr cast 

10 dim x:<Se3 

28 Print chr*<1473 

30 print "OrMtett Common Divisor 

35 print 

40 inPut "Nunbtrs In stf'Jn 

50 Print "Numbtrs"; 

60 s-ltZS 

70 for k-1 to n 

S0 inPut xCk) 

90 If x<k»s thtn 110 

100 s>x(k> 

110 ntxt k 

120 9-0 

130 for n"2 to 5 
140 for i«l to n 

150 if x(iVi>iOint<x(i thtn 180 
les ntxt 1 
170 9«M 
180 next m 

190 Print "NuKibtrs "; 

200 if 9>0 thtn 230 

210 Print "art reUtlvtls Print. 

220 9oto 240 

230 Print "OCD is"; 9 

240 inPut "Another (n.nVStt 

250 if z*-"»" or z«""V" thtn 35 

260 Print "OkiK. But for now." 



Ortattst Conntn Divisor 



WuKibers in set? 3 
Nunbtrs? 12 
36 

Numbers OCD is 4 
Another <«,n)7 k 

Numbtrs in set? 3 
Numbers? 20 

36 

96 

Numbers OCD Is 4 
Another <M.n)? » 

Numbers in set? 3 
Numbers? 20 

36 

97 

Numbers are relativelK Prime. 
Another <s*jn)? n 



Okas*. B«t for now. 
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In the last section we discussed prime numbers. Here is an interesting chal- 
lenge for you involving prime numbers. Until late 1982, the longest pro- 
gression of prime numbers in which all differed by the same number was 17. 
Prof. Paul Pritchard in the computer science department at Cornell Univer- 
sity wrote a program to determine if there was a longer progression. Using a 
DEC VAX-1 1/780, he found the string of 18 numbers shown below. He also 
discovered fourteen other 17-number progressions and ten 18-number pro- 
gressions, but none yet with 19 numbers. He believes there is at least one; can 
you find it? 



107928278317 
117851061187 
127773844057 
137696626927 
147619409797 
157542192667 
167464975537 
177387758407 
187310541277 



197233324147 
207156107017 
217078889887 
227001672757 
236924455627 
246847238497 
256770021367 
266692804237 
276615587107 
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Cryptarithmetic Problems 

Cryptarithmetic or alphametic problems are arithmetic expressions in 
which the digits are replaced by letters of the alphabet. Each digit is asso- 
ciated with a letter to produce an interesting statement, for example: 

SEND 
+ MORE 

MONEY 

If the college student who sent this message to his father needed $106.52 
for plane fare home, this was the right message to send since this combination 
of letters has one unique solution, in particular: 

9567 
+ 1085 

10652 

However, if the student let things go to the last moment and was in more of 
a rush, he might have reworded the message: 

WIRE 
+ MOR E 
MONEY 

In this case, how much should his dad send? Earlier in the book, trial and 
error approaches to solving problems were discussed. It was noted that the 
brute force approach of trying every alternative was sometimes appropriate. 
Is it in this case? 

No! The number of possible alternative solutions is the factorial ot tne 
number of different letters in the alphametic expression, i.e., 8! or 40,320. A 
program to try out every one of these possibilities would run for a 1-o-n-g 
time. 

In this case it is much more efficient to apply some common sense to nar- 
row down the number of aUernatives. The best approach to this process is to 
divide up the search space into large classes (or sets), according to a common 
property shared by members of each class, and then attempt to eliminate en- 
tire classes by the method of contradiction. 

Consider the "WIRE + MORE = MONEY" problem. Can E = 0? Since 
E + E = Y, Y must also equal 0, contradicting the fact that Y and E must be 
different digits. Thus, the entire class of solutions in which E = can be 

ruled out. „ 
Consider E = 3. Now Y = 6 and there is no carry to the next column, bo 

in this column R + R = EorE+10 if a carry is involved. But in either case 

E must be an even number since 2R is always even; this contradicts the 

assumption that E = 3. 

By following this type of classificatory contradiction process for each of the 
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digits in the order E, R, I, O and N, the computer program will search out all 
possible solutions to the problem. Unlike the "send more money" problem, 
the "wire more money" problem has five possible solutions. A smart father 
would choose the lowest solution and wire his son $103.48. 



5 Print chr*<14> :r«m switch tower c»s» 
la Print chr*< 147) 

20 Print "CrKPtarithnstlc »olvtr...." 
30 Print 

40 Print " W I R E" 
S0 Print " + M R E" 

£0 Print " " 

70 Print " « N E V" 
90 Print 
100 m-1 

110 for t«2 to 9 
120 U'fft 

130 If v>10 thtn 150 

140 el-0 

150 Soto 180 

160 cl-1 

170 !*-»-ia 

180 for r«0 to 9 

190 if r-m or r«» or r-a thtn 480 

200 if r+r+cl«t thtn 230 

210 if r+r+cl-t+10 thtn 250 

220 Soto 480 

230 c2>0 

240 Soto 2^0 

250 c2«l 

260 for 1»1 to 9 

270 if i"m or i"t or i«« or i"r thtn 470 
280 for o"0 to 9 

290 if o"iii or o"t or o*v or o"r or o"l thtn 460 

300 n"i+o+c2 

310 if n>«10 thtn 340 

320 c3-0 

330 Soto 360 

340 c3-l 

350 n«n-10 

360 if n"iii or n"t or n"y or n"r or n"l or n"0 thtn 460 
370 for u«0 to 9 

380 if w"n or u"t or or w"r or w"l or w"0 or w"n thtn 450 
390 if o+10<>u+n+c3 thtn 450 
400 Print 

410 Print " ";w;i;r;t 
420 Print " + ";m;o;r;t 

430 Print " » 

440 Print n;o;ri',e!vl 
445 Print :print 
450 next w 
460 next o 
470 next i 
480 next r 
490 next e 
500 »toP 
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CrBPtarlthmefcIc »olwer. 



U I R E 
+ M R E 



M N E V 



9 7 6 2 
+ 1 « 2 



10 8 2 4 



9 5 7 4 
+ 10 7 4 



1 6 4 8 

9 2 8 7 

+ 10 8 7 

1 3 7 4 



9 2 7 4 
+ 10 7 4 



1 3 4 8 



9 3 8 7 
+ 10 8 7 



1 6 7 4 



There are other approaches to solving cryptarithmetic problems, but all of 
them benefit greatly from reducing the search space as much as possible be- 
fore putting the problem on the computer. See if you can devise another 
successful approach and write a program to implement it. 
Here are some problems for you to try. 



DONALD 
+ GERALD 

ROBERT 

ABC 
X DE 



FEC 
DEC 

HGBC 



TWO 
X TWO 

THREE 

THE 
EARTH 
VENUS 
SATURN 
+ URANUS 

NEPTUNE 



ABCDE 
X 4 

EDCBA 

SPRING 
RAINS 
BRING 
+ GREEN 

PLAINS 



ONE 
TWO 
+ FIVE 

EIGHT 



FORTY 
TEN 
+ TEN 

SIXTY 



FIVE 
-FOUR 

ONE 

-I- ONE 

TWO 



VIOLIN + VIOLIN + VIOLA + CELLO = QUARTET 



THREE + NINE = EIGHT + FOUR 
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Sailors and Monkey Problem 



There are many variations of the sailors and monkey problem. Here is one 
of them. 

Five sailors and a monkey were on an island. One evening the sailors 
rounded up all the coconuts they could find and put them in a large pile. Be- 
ing exhausted from working so hard, they decided to wait and divide them up 
equally in the morning. During the night, a sailor awoke and separated the 
nuts into five equal piles, but had one nut left over which he gave to the mon- 
key. He took one pile, hid it, and pushed the other four together and went 
back to sleep. He was followed in this action by the other four sailors, each of 
whom did exactly the same thing. Next morning the remaining nuts were di- 
vided equally with one remaining nut going to the monkey. What is the small- 
est number of coconuts with which they could have begun? 

Although there is an elegant algebraic solution to this problem, a more 
suitable approach for the computer is that of working backwards. A typical 
solution to a problem can be thought of as a path that leads from the given 
information to the goal. However, in this case the goal, or final state, is 
known, thus it is easier to start there and work backwards to the initial state. 

As mentioned at the outset, many sailor and monkey problems exist, in 
fact, an infinite number of them. For example, instead of five sailors, there 
could be three or six or 14. Thus it is desirable to devise a general solution in- 
stead of just one to solve one specific problem. 

In the flowchart and computer program, S is the number of sailors and A is 
the number of coconuts that each sailor received in the final division of the 
pile. Since one coconut was given to the monkey at each division, the total 
number of coconuts left in the morning must be S x A + 1. But this pile came 
from pushing together S - 1 equal piles. Thus, the key condition that must 
hold for (S x A -I- 1) / (S - 1) to be an integer K, which represents the num- 
ber of coconuts that the last sailor stole from a pile of S x K + 1 coconuts. 
But this pile is the result of pushing together S - 1 equal piles by the previous 
thief, so again (S x K + 1) / (S - 1) is an integer and so on back through all S 
raids on the pile. 

Note in the flowchart (and program) that the first trial value for A is 1 
(Statement 40). In Statement 80 this value is increased by 1 until the value of 
(S X K + 1)/ (S - 1) is an integer as tested for in Statement 70. This process is 
then continued until the counter for the second loop, N (nighttime pile di- 
visions) equals the number of sailors. 

Although the program will work lor any number of sailors, it takes a fairly 
long time to run for more than five. Remembering what you have learned ear- 
lier in this chapter, can you devise a way to make the program more efficient? 
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5 Print chrt< 14) :rt« twitch lowtr c«m 

le Print chr«<147> 

20 Print "Stilor* i Monke«" 

30 input "No. »iHor»";t 

40 «-l 

30 k't 

£0 n*0 

70 If <»«k+l)/<i-l).int«»«k+l)/'<»-l)5 then 100 
S0 fti*\ 
90 fioto S0 
100 n«n+l 

110 k"<»«k+l>A»-l) 
120 if n-« th«n 140 

130 9eto 70 

140 Print "Coconutt -";«*k+l 

150 Print "In tht nornlnO, t»ch Mllor 9tt»";«; 



Stilor* i MonkM 
No. ftcilorft? 3 
Coconut! ■ 79 

In th« nominfj t«ch ««llor 9tts 7 



Sailor* i Monkty 
No. Milors? S 
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Super Accuracy 

Under normal circumstances, your computer performs computations to six 
or seven digits of accuracy. Double precision computations increase accuracy 

to 13 digits or so. ■ a ■ 

However, it is possible to do computations one digit at a time and assign 
each digit to an element in an array. This will achieve virtually any desired 
accuracy. Any, up to the maximum array size that is. 

This example program performs the rather simple operation of successively 
doubUng a number (which is the same as raising 2 to a power). 
If the number to be represented is 8192, then: 

A(4) A(3) A(2) A(l) 

8 19 2 

To add this to itself, first the rightmost digits are added: A(l) + A(l). 
If there is a carry, the variable C is set equal to 1, otherwise C is 0. The total 
is put into B in Line 100. 

If B is less than 10, there is no carry (C=0) and the new A(l) equals B. If 
B is greater than 10 there is a carry (C=l) and the new A(l) equals B-10. 

This operation is continued for all the digits (D) of the number and, when 
it is finished, the new number is printed in Lines 190-210. 

If A(N) is printed followed by a semicolon(;) for tight packing, the Basic 
print routine would leave a space in front of each digit (for the sign) and a 
space after each digit (for readability). In the program here, these spaces are 
not wanted, hence the print routine in Line 200 is used which prints the string 
value of the ASCII value of each digit (which is the same as the digit itselO 
but without the spaces. 

This program, incidentally solves the challenge to calculate the number ot 
moves in the Towers of Brahma problem (see "Change for Any Amount to 
$5.00). The approach is also used in the next section, "Palindromes." 
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5 Print chr*<14):r»iii twitch lower case 

10 Print "Computes 2 to Nth to «ccur«cj" : Print 

20 dim t(ie0) 

30 B"0:c"0 

40 d-0 

50 for i-1 to 50:t<i)»0:next i 

60 4l<i)-l 
70 i«0:c-0 
80 ««m+l 
90 i"i+l 

100 b"«<i>+*<i)+c 

110 if b<10 then c"0:9oto 140 

120 b-b-10 

130 c-1 

140 «(l>-b 

iS0 if iCd then 90 

160 if c«l then 90 

170 Print m; 

180 d-i 

190 for n"i to 1 tteP -1 
200 Print chr*<«<n>+4B); 
210 next n 
220 Print 
230 Ooto 70 



ConPutes 2 to Nth to tny accuracy 



1 


2 


34 


17179869184 


2 


4 


35 


3435973S368 


3 


8 


36 


68719476736 


4 


16 


37 


137438953472 


5 


32 


38 


274877906944 


6 


64 


39 


S497SS8 13888 


7 


128 


40 


1099511627776 


8 


2S6 


41 


2199023255552 


9 


512 


42 


4398046511104 


10 


1024 


43 


8796093022208 


11 


2048 


44 


17592186044416 


12 


4096 


45 


35184372088832 


13 


8192 


46 


70368744177664 


14 


16384 


47 


140737488355328 


IS 


32768 


48 


281474976710656 


16 


65536 


49 


562949953421312 


17 


131072 


50 


1 125899906842624 


18 


262144 


51 


2251799813685248 


19 


524288 


52 


4503599627370496 


20 


1048576 


S3 


9007 1 99254740992 


21 


2097152 


54 


1 80 1 439850948 1 984 


22 


4194304 


55 


360287970 1B963968 


23 


838B608 


56 


72057594037927936 


24 


16777216 


57 


1441151 88075855872 


2S 


33554432 


58 


28823037615171 1744 


26 


67108864 


39 


576460752303423488 


27 


134217728 


60 


1 152921504606846976 


28 


268435456 


61 


23O58430092136939S2 


29 


536870912 


62 


46 1 1 6860 1 8427387904 


30 


1073741824 


63 


9223372036B54775808 


31 


2147483648 


64 


1 844674407370955 1616 


32 


4294967296 


65 


36893488 147419103232 


33 


85B9934592 


66 


73786976294838206464 
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Palindromes 

A palindrome is a word, verse, or number that reads the saine backwards 
ortorwal For example, the words "mom" and "eye" are palmdromes. So 



are each of the lines in this verse: 



Egad, a base life defiles a bad age 
Doom an evil deed, liven a mood 
Harass sensuousness, Sarah 
Golf; No, sir, prefer prison-flog 
Ban campus motto, "Bottoms up, MacNab 



Numeric palindromes are those numbers which read the -me^^^^^^^ 
forward. The examination of these numbers is a field rich with possibilities 

^^^5::::™^^^^ pal— r^ses an — — 
question. Begin with any positive integer. If it is not Y^-j^T. 
digits and add the two numbers. If the sum is not a Pf "drome- ^^^^^^^^^^ as the 
original number and continue. The process stops when a palindrome is 
Vnr pxnmnle. beeinning with 78: 



uiigiiiai iiuiiii^vi "-iv ~ 

tained. For example, beginning with 78 

78 
+ 87 

165 
+ 561 



726 
+ 627 

1353 
+ 3531 



The conjecture often assumed true, is that this process will always lead to 
a patdX Arid indeed that is Just what usually happen. Mo«be. 
less than 10,000 will produce a palindrome m l^^^f^^Jl^^^ 
there is a real thorn in the side of this conjecture, the number 196. ^^n you 
determineTf a palindrome will ever be produced with a starting number of 

1675? Because ten other "umbers under 1000 wiU^^^^^^ 
1675 and thus may not become palmdromic. The tirst live oi 
are 196. 295, 394. 493. and 592. What are the other five.' 
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The program here will accept any number as a starting value and complete 
the process of adding the successive reversals and testing if the sum is a palin- 
drome. Try it with some numbers and see if you can identify any patterns. 



5 Print chrt<H):rti<i iwltch lowtr c«tt 

10 Print "TMt» conJtcturt thit uddlnfl ("•vtritd di8lt« will Product « Palindron* 
20 Print 
30 dim b<Se> 
40 Print 

50 InPut "Enttr numbtr";* 
€0 t-0 
70 t-t+1 
80 t-i/'ie 

90 If lnt<t»0 th»n 70 
100 for c"t to 1 »ttP -1 

110 «>««10 

120 b(c)«int(«-10itilnt<«/'10)} 
130 ntxt c 
140 d-0 

150 for c-1 to int<t/2> 
1S0 If bCc)-b<t+l-c> then 180 
170 d"l 
180 ntxt c 

190 for c"t to 1 »ttP -1 

200 Print chr«<b<c>+48>; 

210 ntxt c 

220 if d-1 thtn 260 

230 Print " P*Undro»tl " 

250 «oto 40 

2£0 Print " Not utt" 

270 if t/'2<-lnt<t/'2> thtn 290 

280 b<lnt<t^2)+l)"2Kib<lntCt/2>+l) 

290 for c«l to intli/Zy 

300 b<c>-b<c>+b<t+l-c) 

310 ntxt c 

320 for c"l to lnt<t^'2> 

330 b<t+l-c)«b<c> 

340 ntxt c 

350 bit+D'H 

360 for c"l to t 

370 b(c+l)-b<c+l>+lnt<b<c)<'10) 

380 b<c)"b<c>-10«lnt(b<c>/'10) 

390 ntxt c 

400 if b(t+lK>0 thtn 140 
410 t-t+1 

420 goto 140 



Tt«t» conJtcture that «ddlnfl rtvtritd di91t» will Product * Pulindront 



Enttr Numbtr 19 
19 Not Ktt 
110 Not ytt 
121 Pallndront! 
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Enttr Number 96 
96 Not ytt 
1£S Nofc ytt 
726 Not Ktt 
1353 Not xtt 
4884 PtUndroMt! 



Enttr Number 196 
196 Not vet 
887 Not yet 
1673 Not yet 
7436 Not yet 
13783 Not yet 
32514 Not yet 
94839 Not yet 
187088 Not yet 
1067869 Not yet 
10753470 Not yet 
18211171 Not yet 



Using this method, write a program that exammes all integers bet^^^^^^^^ 
1 and 10.000 excluding those that sum to 1675 at 7 P^-^^Jf^^^^^^^^^^^^^ 
show? By the way, you will have to devise a way to deal with 14-dig.t integers 
which are larger than your computer can normally handle. 

Huh? Is this program "too hot to hoot?" 
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4 

Convergence and Recursion 



The computer is especially suitable for doing repetitive and tedious calcula- 
tions. Two mathematical approaches for solving problems that involve repet- 
itive calculations are convergence and recursion. 

Some problems can be reasonably easily stated in words or described with a 
few simple equations but there are many possible solutions. For example, how 
many ways can you make change for a dime? It is simply stated and the num- 
ber of ways can be enumerated fairly easily: two nickels, one nickel and five 
pennies, or ten pennies, three ways in all. But if you want to solve for all the 
ways of making change for a dollar or five dollars, it would be nice to have 
some help. 

Help on this kind of problem comes from a class of computer program that 
simply breaks the problem into smaller ones and counts up all the alternative 
solutions according to a set of rules. But an even more powerful technique is 
known as recursion. Using this technique, a simple solving algorithm or rou- 
tine is set up to solve the smallest subset of the problem. The unique power in 
a recursive routine comes from the ability of the routine actualUy being able 
to call itself. This is discussed further in the second program in this section. 

Another approach for solving problems that do not have an exact answer is 
that of successive approximations. For example, the exact value of pi, e or the 
length of an irregular curve cannot be precisely determined. But by means of 
increasingly accurate approximations, it is possible to approach the desired 
value from above or below or to converge on it from two directions. The last 
four programs in this chapter illustrate successive approximations and 
convergence. 
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Change For a Dollar 

Even though there is very httle you can buy for a penny these days, the 
coin will probably be around for some time to come since it is needed to make 
change for odd amounts of sales tax and to fill up penny collections. 

Today U.S. coinage consists of five coins: penny, nickel, dime, quarter, and 
half dollar. How many ways can coins of these denominations be used to 
make change for one dollar? For example, one way is two half dollars, an- 
other is one half dollar and two quarters, and so on. Make a best guess now 
and write it down before you read further. 

There are several different ways to approach a problem of this kind. One is 
to break it down into smaller, nore easily solved problems. In other words, 
how many ways can you make change for a quarter? For a dime? You would 
solve these subproblems and combine the answers to give the overall solution. 

If you were more mathematically inclined, you could write a series of equa- 
tions relating each piece of change to every other one and to the dollar and 

solve them. , . • t i 

A third approach is to do the problem by writing down combinations until 
all the different possibilities are exhausted (or until you are exhausted) and 
then count them all up. This might be called solving the problem by exhaus- 
tion and is a method quite suitable for putting on the computer. 

Write a program that uses this approach to solve the problem. If you use 
loops and count by one, it could take a long time for the computer to run 
through all the possible combinations, possibly many hours. 

Also, if you want to print out all the possible combinations, be warned that 
the printing could also take quite some time and a fair amount of paper. 
There are more combinations than you might think! 

In fact, most people will not be able to guess the answer to this problem, or 
even come close. Ask several of your friends how many ways they think a 
dollar can be changed. Record all the responses and then tabulate them on 
your computer. What is the mean (average) of all the guesses? The extremes? 

The program included here uses the first method discussed to solve the 
problem, in particular, breaking down the problem into subproblems and 
then combining the solutions into one final answer. 

First, the main problem is broken into the next smaller one of making 
change for half dollars. There are three such problems: no half dollars (H - 
0), one half dollar (H = 1), and two half dollars (H = 2). The last problem is 
trivial since there is only one way, but the other two need to be broken down 
further. 

This is done by dividing the remaining money into quarters and consid- 
ering the subproblems on down to the lower denominations. As the number 
of subproblems is expanded, each one becomes easier to solve. In fact, 
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subgoals, which can be solved in only one way, are finally reached. For exam- 
ple, if H = 1, Q = 1, D = 2, and N = 0, then the pennies (P) must equal 
five in order that the total equal 100. 

Notice that at the quarter, dime and nickel stages, adjustments are made in 
the limits of the loops depending upon how much money there is left to 
change. For example, if H = 1, the only possible subgoals for quarters are 0, 
1, and 2, but not 3 or 4. Also notice that there is no need to test combinations 
of coins to see if they add up to 100, nor is it necessary to include the penny 
as a variable. Simply counting the number of subgoals is sufficient since each 
one can be solved in only one way. 



3 Print chr*<14) :rt(« twitch lower c*»» 
10 erne 

20 for h"0 to 2 

30 for <1"0 to 4-2ilih 

40 for cl-0 to l0-5llih-2.5W 

50 for n-0 to 20-10«h-5#q-2«d 

60 c-c+1 

70 next n 

80 next d 

90 next 4 

100 next h 

110 Print "fl dolUr c»n bt chanfled in";c; "different w*»»." 



H dollmr c«n be ch*n9ed in 292 different uasit. 



Try to make some changes in this program or write a new one to solve the 
following problems. Say you want two quarters in your change to play some 
video games. In how many ways can a dollar be changed to provide at least 
two quarters? 

Visiting a small town, you find the parking meters still take pennies. In 
how many ways might you get change so that you had at least three pennies? 
Is this any different than the number of ways that would give you five pen- 
nies? Say you want to make a phone call also; in how many ways can you 
change a dollar to produce at least four pennies and one dime? 
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Change For Any Amount to $5.00 

Another way to attack the change problem in the previous section is by 
means of the programming technique called recursion. Get famihar with this 
one-it is very powerful! Donald Piele and Larry Wood described this 
method in an issue of Creative Computing. 

First, define the variables which represent the number of ways to make 
change for n cents using the coins specified; 
A Only pennies 
B Nickels and pennies 
C Dimes, nickels, and pennies 
D Quarters, dimes, nickels, and pennies 
E Halves, quarters, dimes, nickels, and pennies. 
Initially, there are two subproblems in making change for n cents. In the 
first, no half dollars are used, and D is the number of ways to change n cents. 
Second, when one or more halves are used, after one is paid, there remain 
n-50 cents to pay which can be done in En-50 ways. . p _ 

Since these two cases are mutually exclusive, it can be inferred that E„ - 
Dn + E„_5o. Similarly, 

D„ = Cn + D„_25 
Cn = Bn + Cn-10 

Now, begin with" the simplest case and build up to E,oo. First o^all, it is 
easy to understand why E^ = 1. From above, when n = 50, E.^ - D50 + 
E„, and it is possible to make change for 50 cents only one more way if^alt 
dollars are allowed. Therefore = 1. Likewise, D^ = = -A^ - l- 
It is also true that An = 1 for all values of n since there is only one way to 
make change using only pennies. Now the recursive relationships can be used 
to solve the original problem. 

This is the strategy used in the program. It also has the added advantage 
that it can count the number of ways of making change (with corns) for any 
specified amount. 

Now it is your turn. Can you modify the program here to include dollar 
bills so it could count the number of ways to make change for any amount up 
to $10.00? . 

Using any method of change making you prefer, write a program to make 
change for one ruble. Russian coins come in denominations of 1, 2, 3, 5, 10, 
15, 20, 50, and 100 kopecks. There are 100 kopecks in one ruble. 
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5 Print chr«<14>:rtn twitch loucr cut 

18 din ii<iei>;b(iei)^c(iei)<di;iei>>*ciei) 

13 Prlflt 

20 InPut "flmount to bt ch«n9td";x 
30 n-lnt<2mx)+l 

40 *Cl>«i:b<l)-l!c<l>"i:d<l>-l!t<15"l 
50 for J"2 to n 
£0 «CJ>>1 

70 b<J)"*<J)+b<J-l) 

80 c(J>«b(J) 

90 If JC-2 thn) 110 

100 c<J>«b<J>+c<J-2> 

110 d(J>c(J> 

120 If J<-5 th*n 140 

130 d<J)"c<J)+d<J-5) 

140 ((JJ-dCJ) 

150 If K-10 thtn 170 

160 t<J>-d<J>+t<J-10> 

170 ntxt J 

180 Print "Vou c*n m«kt ch«nM lin";f<m);"u«Dt." 
190 goto 15 



Rmount to bt ch«n9td? 1.13 

Vou can makt ch«n«t in 384 u«y«. 

flmount to bt chinfltd? 4.25 

Vou ctn mtkt chtnSt In 32985 u«ys. 

Amount to bt chunStd? .1 

Vou c«n makt ch»n9t in 4 iii«»s. 

Rmount to bt ch«n9td? .12 
Vou can makt chanSt in 4 Mm. 
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Perhaps the most famous problem used to demonstrate the pnnciples of 
recursion is the Towers of Brahma. It is sometimes called the Towers of Ha- 
noi or Pharoah's Needles. Here is the problem in as close to ongmal form as 
possible. You should be able to solve it with a relatively short program usmg 



recursion. 



In the great temple at Benares beneath the dome which marks the center ot 
the world rests a brass plate in which are fixed three diamond needles, each a 
cubit high and as thick as the body of a bee. On one of these needles, at the 
Creation, God placed 64 discs of pure gold, the largest disc restmg on the 
brass plate and the others getting smaller and smaller up to the top one. This 
is the Tower of Brahma. 

Day and night unceasingly, the priests transfer the discs from one needle to 
another, according to the fixed and immutable laws of Brahma. These laws 
require that the priest on duty must not move more than one disc at a time 
and that he must place this disc on a needle so there is no smaller disc below 
it When the 64 discs shall have been thus transferred from the needle which, 
at the Creation, God placed them, to one of the other needles; tower, temple, 
and Brahmans alike will crumble into dust, and with a thunderclap, the 
world will vanish. 

If the priests were to effect one transfer every second, and work 24 hours 
per day for each day of the year, it would take them 58,454,204,609 decades 
plus slightly more than six years to perform the feat, assuming they never 
made a mistake— for one small slip would undo all the work. 

How many transfers are required to fulfill the prophecy? Try out your pro- 
gram with fewer discs than 64 to make sure you are on the right track. Here 
is a table of the first few transfers: 



Discs 


Moves 


Discs 


Moves 


1 


1 


6 


63 


2 


3 


7 


127 


3 


7 


8 


255 


4 


15 


9 


511 


5 


31 


10 


1023 
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Converge on e and Pi 



An incredibly important mathematical constant is designated by the small 
letter e. This constant is both irrational and transcendental. Look up those 
terms in a dictionary or math book if you wish, or just plunge on to the next 
paragraph. 

The constant e was first derived by John Napier, also the inventor of loga- 
rithms, to whom we owe an eternal debt of gratitude. Why? If e had never 
been discovered, advances in mathematics, physics, and astronomy would 
have lagged a century or more, because e is the base of all natural logarithms 
and these logarithms are the basis for many branches of science and 
mathematics. 

How is e calculated? The constant e is the limiting value of this expression 
as n approaches infinity: 

e = (l + yn)n 

Its exact value can never be found, but to 15 places e equals 
2.718281828459045 . . . How can e be calculated? First take l-Vz and square 
it; that equals 2-74. Then cube l-'/g and you get 2.3686. Raising l-Vj to the 
fourth power gives 2.414, and so on. Write a program for this method and 
have it print out the initial value of e and each value after each additional 
fraction is added on. 

Another approach is to expand the expression above using the binominal 
theorem and, again, letting n approach infinity. The expression for the expan- 
sion is: 

e=l + l + 2! + 3l"*'4l + r!---;?T 

Now here is where the computer can again be of some assistance. Since 3! is 
3*2! and 4! is 4*3!, all the calculations need not be done for each additional 
fraction. Look at the program and particularly note the calculations in State- 
ments 70 and 90. 



5 Crlnt chr«<14}:rtm twitch lower cum 

10 Print "Converse on e" 

20 e-1 

30 1-0 

40 1-1+1 

50 d-l 

50 forJ-1 to 1 
70 d-d«J 
80 next J 
90 e-e+l/'d 
100 Print e 
110 »eto 40 



59 



Converse on e 

2 

2.5 

2.6666£fiS7 

2.70833334 

2.7166fi£S7 

2.71805556 

2.71825397 

2.71827877 

2.71828153 

2.7182S18 

2.71828183 

2.71828183 

2.71828183 

2.71828183 

2.71828183 



Pi is another important mathematical constant that is irrational (meanmg 
its exact value can never be determined) and transcendental (meaning it is not 
the solution to any algebraic equation). Interestingly, pi was known and used 
by the ancients. Archimedes, who lived in the second century B.C., by using a 
regular polygon of 96 sides (nearly a circle), proved that the value of pi was 
less than 22/7 and greater than 3-^%i, a remarkable achievement for the 
mathematics of his day. 

Ptolemy in 150 A.D. used the value of 3.1416 for pi and in the middle ot 
the sixteenth century the amazing fraction 355/113 was discovered, giving 
the value of pi accurately to six decimal places. 

Incidentally, in 1897, the General Assembly of Indiana passed a bill ruhng 
that the value of pi was four. 

Several infinite series can be used to grind out increasingly accurate values 
for pi. One such series is (l-Va + Vs - Vt + '/b • • • )• This series is called an 
arithmetic series and converges very slowly. The program here displays only 
every 500th value of the series. Don't be alarmed if the program does not 
seem to be running very fast; a fair amount of calculating is going on between 
each value that is printed. 



5 Print chr«<14):reKi switch lower c»»* 

10 Print chrf<147) „ 

20 Print "Conv»r»t on Pi flrlthmetic »trie» / 

30 t-1 

40 1-1 

50 

60 P«0 
70 q-fl+1 
80 P"P+t/'i 
90 i"l+2 
100 »"-» 

110 if fl<499 thtn 70 
120 1-0 

130 Print:Print P#4; 
140 eoto 70 



Conwers* on Pi firithmetlc strle* 
3.14359667 
3.14059066 
3. 14226067 
3.14109167 
3.14199346 
3. 14125867 
3.14187896 
3.14134218 
3.14181534 
3.14139227 
3.14177485 
3. 14142568 



60 



Actually, the approach used by Archimedes converges much more quickly 
and, with the aid of a computer, it is possible to go far beyond the 96-sided 
polygon used by Archimedes. 



His approach was to construct inscribed and circumscribed polygons and 
measure the perimeters to approximate the circumference of a circle. Con- 
sider a polygon circumscribed around a circle of radius 1. The perimeter 
equals the length of one side times the number of sides. Since the tangent of x 
= AB/BC, but BC = 1, then tan (x) = AB and the length of a side = 2 tan 
(x). Since the circumference = 2 pi r and r = 1, then pi is the circumference 
(or perimeter of an n-sided polygon) divided by 2. 

Similar trigonometry leads to the perimeter of an inscribed polygon being 
equal to the number of sides times sin (x) * cos (x). 

The second program produces values for pi using inscribed and circum- 
scribed polygons. Unfortunately, there is one large flaw in the program, be- 
cause degrees must be converted into radians in Statement 40. This means, of 
course, that you must already know the value of pi, since the conversion fac- 
tor is 360 degrees divided by 2 pi. 

Setting this flaw aside, it is interesting to note how quickly this program 
converges on the value of pi compared to the preceding one. That is because 
this one converges geometrically rather than arithmetically. 

Can you figure out a geometric convergence to the value of pi that does not 
require that you know it (or a conversion factor) before you start? 

S Print chrt<14);r«iii switch lower case 
10 Print "Converse on Pi ba PolaSons. " :prlnt 
15 Print" Inscribed Circumscribed" : Print 
20 n-£ 
30 n«2Kin 




40 x"3S0^<n«57. 29578) 
50 Print n«»in<xJ«cos(x)/'2; 
S0 Print t«b<13);n«tan<x)^2 
70 9oto 30 



Converse on Pi bs PotsSons. 



Inscribed 



Circumscribed 



2.5980762 
2.99999998 
3. 10582851 
3. 13262859 
3.13935018 
3.14103193 
3. 14145245 



3.46410158 
3.21539028 
3. 15965992 
3.14608S19 
3. 14271458 
3. 14187303 
3. 14166272 
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Convergence on Pi Revisited 



In answer to the question posed in the last paragraph of the preceding sec- 
tion, here is a way to converge on pi without knowing its value beforehand. 

As in the previous program, the basic approach is to add up the length of 
the sides on an inscribed polygon and divide by 2r to obtain a value for pi. 
The program starts with a square (four sides) and doubles the number of sides 
each time. If the old side length is S, then the length of S' of a side of a new 
polygon with twice as many sides is obtained by applying the Pythagorean 
theorem. In particular, 

X^' + (S/2)^ = 

(R-Xf + (S/iy = (S')' 

Thus, J — — 

S' = V(R-VR^-(S/2)T + (S/2f 
It is easy to reduce this formula algebraically, but accuracy suffers if this is 
done. Also, you will find that S*S is slightly more accurate than St2. 

Unfortunately, striving for maximum accuracy is somewhat moot on a 
computer that does not have double precision arithmetic. Notice that ac- 
curacy does not improve with more than 1024 sides and, indeed, as numbers 
in the calculations start to exceed the capacity of the computer (4194304 
sides), the accuracy starts to deteriorate badly. 

There is yet another method to compute pi by convergence. It uses discov- 
eries of Gregory and Euler. Gregory discovered the formula for arctangent: 
Arctan x = x - xV3 + xV5 - xV7 + ... 
Euler came up with a rather interesting formula for pi: 

pi = 4 (arctan CA ) + arctan (Vis)) 
See if you can combine these two formulas to calculate pi. If you are very 
clever, you can do your calculation to yield far more than the seven decimal 
place accuracy obtained by the programs presented here so far. 




5 Print chr«<14>:rtM lultch lower cue 

10 Prlflt "Convtrse on Pi bu Infcrlbed l>oly9on»" 

20 Print "Sidtt Ptrimtttr" 

30 r-ie 

40 z-2 

50 n"S+z 

60 »"rlll»qr(i> 

70 for k"l to 28 

80 Print n;t«ba0);<t«n>/'<z*r) 

90 ii«sili«/'<zi|iz> 

100 x-r-%<iri:r»r-v) 

110 »"«<Ir<x«x+») 

120 nartiDz 

130 ntxt k 



Convfrs* on Pi bu Inscribtd Pol»9ont 



Sidti 


Perlrntttr 


4 


2.82842713 


8 


3.06146746 


16 


3.12144515 


32 


3.13654849 


64 


3.14033116 


128 


3. 14127725 


256 


3.1415138 


512 


3. 14157294 


1024 


3.14158773 


2048 


3.14159142 


4096 


3. 14159235 


9192 


3. 14159258 


16384 


3.14159264 
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Length of Any Curve 



The previous programs have demonstrated how it is possible to compute a 
very accurate value of pi by adding together the length of the sides of a poly- 
gon as it approaches a circle and dividing by 2r. 

Using a similar approach, it should be possible to inscribe a polygon, or 
portions of a polygon, inside any regular curve and thus determine the length 
of the curve. This program approximates the length of any curve as defined in 
Statement 100 by dividing it into an increasing number of subintervals and 
computing the sum of the secants (a straight line that cuts a curve at two or 
more points). 

To run the program, you must enter the formula or equation describing 
your curve in Statement 100 in the form: 

100 DEF FNA(X) = your function of X 

You then type RUN and enter the end points of the curve you want to use 
in your calculation. These points are entered in the form of the abscissa, 
which means the horizontal (or y) coordinate of the point. 

The program is written to sum the successive secant lengths and to cal- 
culate the percent of change in each summation compared to the preceding 
one. The sample run uses the function 2x^ + 3x^ -2x +3. Note the substan- 
tial improvements in the length calculation as the number of intervals in- 
creases from 2 to 16, but the rather slight improvements beyond 16. 

Try this program with different curves and functions. It might help to plot 
the function first (remember the program to do that?) and then compute its 
length. Is this method of length calculation more accurate for a function with 
no changes in direction of the curve within the interval selected or for a func- 
tion with one or more changes? Why? 
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3 Print chr*<14):rem switch lower ca&e 

10 Print chr*<147) 

20 Print "Curve length." 

30 Print "Define function in line 100" 

40 Print 

50 Print 

60 inPut "flbicis*** of end Polnts"JPifl 
70 Print 

80 Print "Inter Secant ch«n9e" 
85 Print "V«ls LenSth in len9th" 
90 Print 
95 sl>0 

100 def fn*<x>-2«xt9+3«x12-2#x+3 
110 for n«l to 9 
120 •«2t<n-l> 
130 h-<S-P)/e 
140 s>0 

150 for i«0 to e-1 

160 l"»qr<<fn«<P+l«h+h>-fn«(P+l#h))t2+h«h) 
170 s»$+l 
180 next i 

190 if sl>0 then 220 

200 Print e;t*b<5>;»;" No Previous value" 
210 goto 240 

220 P5><<abs<3l-»)>/'ftl>iliie0 
230 Print e;tab<5);»;tab<75;P5 
240 ftl>s 
250 next n 



Curve lenflth. 

Define function in line 100 



fib»ci«sa» of end Points? -1 j 6 



Inter 


Secant 


H chanse 


Valt 


Len9th 


in lenath 


1 


525.046666 


No Previous value 


2 


525. 158263 


.0212546484 


4 


529.652248 


.855739211 


8 


531.017134 


.257694752 


16 


531.964256 


. 178360033 


32 


532.016566 


9.83336555e-03 


64 


532.041679 


4.72022717e-03 


128 


532.048562 


1.29376759e-03 
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Converge on a Square Root 

Since most square roots are irrational, methods used to calculate th^m usu- 
ally involve successive approximations. Although you can simply call up the 
square root function in Basic or on many pocket calculators, it is interesting 
to explore various methods of calculating square roots without these built-in 
functions. After all, these built-in functions are nothing more than successive 
approximation routines already installed in the machme. 

Ooviously, a square root is the inverse of the operation of squaring a num- 
ber All of the methods of calculating square roots use this fact, but the way 
in which it is employed is quite different in various calculators and 

'°The program here calculates an upper and lower limit for the square root 
of a number and, by successive approximations, pinches the root to withm a 
smaller and smaller interval until it reaches the desired level of accuracy. 

The starting value for the lower limit is and for the upper limit the num- 
ber, Z, whose square root is sought. The program then divides this interval 
into ten steps by simply dividing the difference between the nunibers by 0. 
The variable I is increased from the lower limit to the upper one by the value 
of the step, S. At any point, if I squared becomes greater than Z, a new upper 
limit is set to I and a new lower limit is set to I - S. , • , 

This method converges very quickly and adds approximately one decimal 
place of accuracy with each pass beyond the third. What happens when you 
enter into the program a number that has an exact square root such as 25 or 
49? Why? 

Another approach to calculating square roots by successive approxima- 
tions is to start with a trial root, X. If X * X is less than the original number 
N. then increase the trial value by a 0.1. If X * X is greater than N, return to 
the previous value. This is the first digit of the root. Now, start advancing by 
0.01. Continuing in this way, one digit is developed at a time until the desired 

precision is reached. . 

This method is quite suitable and fast for square roots of numbers less than 
1 A good first trial root value is 0.1. But is it suitable for larger numbers? 
How should a starting trial value be determined? In this method, especially 
for numbers greater than 10, the initial trial value for the root matters a great 
deal in determining the length of time it will take for the calculation to con- 
verge. Write a program using this method and compare the speed and ac- 
curacy with the program in the book. 
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3 Print chr«<14):reiii switch lowtr c*»t 
la Print "S<IuJir* Root*" 
20 inPut "Vour nunbtr";z 

40 Print "Lou Limit UP Limit" 
50 «"0 
S0 b-i 

70 »«<b-«)/'10 

80 Print «;t«ba0);b 

90 if »bs<*«b-z><« thtn 180 

100 for i"* to b »teP « 

110 If z<l«l thtn 150 

120 ntxt i 

130 b>b*10 

140 Soto 70 

150 b=i 

160 »■!-« 

170 Soto 70 

180 Print "flccur»t« to .00001" 
190 Print "flvtr«9e "";<»+b)^2 



S<lu*re Root* 
Vour number? 54 
Low Limit UP Limit 

54 

5.4 10.8 

7.02i 

7.344 7.398 
7.344 7.3494 
7.34832 7.34886 
7.348428 7.348482 
7.3484658 7.3484712 
7.34846904 7.34846958 

flccurite to .00001 
flv»r*9e - 7.34846931 



S<lu»r9 Roots 

Vour number? 3 

Low Limit UP Limit 



l./l i.f 
1.731 1-734 
1.7319 1.7322 

flccurj.te to .80001 
flv«ra9« " 1.73203 
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5 

Compounding 



As with successive approximations and recursion, compounding requires 
many repetitive calculations. Some compound interest and growth situations 
can be represented by a formula, but for many problems, solving by repetitive 
calculations is an excellent approach. 

There is nothing magical about compounding; you generally start with an 
initial amount of money, number of animals, etc. This quantity then grows or 
diminishes by a certain percentage at set intervals. The new amount is the old 
amount increased or decreased by this percentage. Repeat this calculation 
over and over again, and you have solved the problem. 

Five different types of problems which involve compounding of some sort 
are described in this chapter. Try the extra problems that follow some of the 
programs; you may be surprised at the results. 
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Indians and Interest 



Here is a simple compound interest problem that produces an astonishing 
answer. The problem has to do with the sale of Manhattan Island to the 
Dutch for approximately $24 worth of trinkets and beads. 

If the $24 that the Indians received in 1626 had been deposited in a bank 
paying 5-y4% interest compounded annually, how much would it amount to 



The program here solves this little problem by making use of the formula 
to calculate compound interest. In particular, if P dollars are invested at an 
interest rate of R (expressed as a decimal) and compounded N times, then the 
total amount A is given by the formula: 



How much was gained in 1983 alone? How much was gained in the decade 
from 1974 to 1983? You can change the value of N in Statement 30 to get the 
answers to these questions. But a better way might be to enter the ending year 
with an input statement. 

Can you read a number expressed in the E (exponential format)? In a more 
conventional format, the number is $11,176,500,000 or $11.2 billion. 

There are many other ways you can improve the program and make it 
more suitable for general purpose compound interest calculations. Modify it 
to accept any starting and ending year, any rate of interest, and any starting 
principal amount. 

The problem as stated is somewhat unrealistic since banks were not paying 
interest rates of 5-y4% from 1626 to 1983. Change the program to calculate 
the total amount based on the following interest rates: 



in 1983? 



A = P(l + R)N 



Years 



Interest Rate 



1626-1830 
1831-1870 
1871-1910 
1911-1921 
1922-1929 
1930-1940 
1941-1945 
1946-1960 
1961-1980 
1981-1983 



1.5% 
2.0% 
3.0% 
3.5% 
6.5% 
2.3% 
3.5% 
5.3% 
6.5% 
9.5% 
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10 P-24 
20 r-.0575 
30 in>19S3-l£26 
40 t«P«<<l+r)tn) 

50 Print "Indisn* Sot «24 In 1626 for NV." 
60 Print "Bk 1983 «t 3.73X it woutd b» *";«; 



Indl»n» 9ot *24 in 162S for NV. 

B» 1983 «t it would b* * 1. 11764926*-i'10 




Systematic Savings 



In the previous program, compound interest was calculated by means of a 
formula well known to bankers, money lenders and real estate agents. How- 
ever, if you did not know the formula, how would you calculate interest on 
money in a savings account by hand or with a calculator? 

You would probably begin by multiplying the principal amount P by the 
interest rate R and adding that amount to the original principal at the start of 
the second year. Doing that for all the years the money is in the bank will 
yield a final amount. Why not write a program to perform the calculations in 
this manner rather than use a formula? Here is such a program. 

It is set up for an initial principal amount of 100 (Statement 50), an interest 
rate of 10% (R = 0. 1 in Statement 60) and ten years (N = 10 in Statement 
70). Naturally these values could be read in using input statements. The 
clever calculation in Statement 100 rounds off the amount to two decimal 
places (dollars and cents). 

In contrast with the program in the previous section, this one does not use 
a compound interest formula, but simply adds the interest each year to the 
growing principal amount in Statement 90. 



5 Print chr«<14>:rtiii twitch lowtr c«»t 

18 Print chr«<147>;"C«lcul«tM Inttrest on *100 Invtsttd lit 10!i for 10 sit»rt" 



20 Print 
30 Print 
40 Print 
50 P-100 
60 r-.l 
70 n-10 
S0 for l«l 
90 P"P+P«r 
100 Print i;t<b<i0). 
110 nixt I 



"fit tnd 
"of »e«r 



to n 



Currtnt" 
b«Unci" 



lnt<P*100*. 5)^100 



CtlcuUtt* Interest on *100 
Invested at IZ'/i for 10 s»e«rs 



fit end 
of i/ear 

1 

2 

3 

4 

5 

6 

7 

S 

9 

10 



Current 
b<Unce 
110 
121 
133.1 
14S.4i 
161.03 
177.16 
194.87 
214.36 
233.79 
239.37 



Now, how could this program be modified to allow for a plan of systematic 
saving? In other words, instead of letting the $100 lie around all lonely while 
it is compounding, each year you add another $100 to it. With this program, 
making the modification for systematic savings is easy: Statement 105 is 
added to add the new deposit each year to the ever growing principal. 
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3 Print chr*<14>:rti» twitch lowtr c«»e 

10 Print ehr*(i47); "C«tcul«te« interest on Invested at every near"; 

20 Print "for 10 »e«r«." 

30 Print "End Current" 

40 Print "Vr Invested B*l» 

90 P-100 

59 c-iee 

60 r-.l 
70 n>10 

60 for i«l to n 
Se p-P+PHir 

100 Print l;t»b<4);iil(c;tib<in;lnt<P«100+.3>^l00 

109 P-P+c 

110 next i 



Calculates Interest on 4100 invested at 10?i every vearfor 10 years. 



End 




Current 


Vr 


Invested 


Bal 


1 


100 


110 


2 


200 


231 


3 


300 


364.1 


4 


400 


910.S1 


9 


900 


671.36 


6 


600 


848.72 


7 


700 


1043.S9 


8 


800 


1257.95 


9 


900 


1493.74 


10 


1000 


1753. 12 



In an effort to attract depositors, many banks over the last 20 years have 
started offering savings and investment accounts in which the deposits are 
compounded more often than annually. In the early 50's, many banks went to 
quarterly compounding; in the late 50's, to daily; and in the early 60's, some 
"competitive" S&L's went to continuous compounding. 

How much does more frequent compounding really mean? Modify either 
of the two programs to compute the interest for more frequent compounding. 
What is the difference in interest for one year for annual, quarterly, monthly, 
daily, and continuous (every second) compounding on a principal amount of 
$1000 invested at 8%? How about for ten years? 

Incidentally, if you want to use the formula method in the previous section 
for this calculation, the formula for P principal invested at R rate com- 
pounded N times per year is: 

A = P(l + R/N)N 

Try this problem which uses the same principles of compounding. Con- 
sumer prices rose an average of 8.8% in 1980. While the government keeps 
trying to bring inflation under control, they don't seem to be meeting with 
much success. Assuming that prices continue to go up this much (8.8%) ev- 
ery year, how much will a $6000 economy car (1980 dollars) cost in the year 
2000? How much will it cost when you are 65 years old? 
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Systematic Savings Revisited 



Using the formula for compound interest and systematic savings, this pro- 
gram will calculate the amount accumulated after a given period of time. 

When you run the program, it will ask how much you wish to save each 
month, the number of compounding periods in a year, the interest rate, and 
the length of time you wish to continue your systematic savings program. The 
program will then calculate the total amount at the end of that period. 

From the preceding programs, you should be able to see that a systematic 
savings program is a very effective way to accumulate a nest egg for the 
future. 



5 Print chr*<14>:rt« twitch lower CMt 

10 Print "C«IcuI«ttt lTtttrt»t and btUnct for t Mtt«m«tic ttvinss Pro9r«n." 
13 Print 

20 input "Monthly dePotlt «";« 

40 InPut "Compounding Ptrlods Per year";b 

80 InPut "Interest r«te <>o<.x>";r 

100 rmr/m 

110 inPut "No. or ve«r<";n 
130 Print 

140 Print "End Current" 

150 Print "Ve»r Inve»ted B*l«nce" 

160 Print 

170 c-12«t 

180 P>12l«t 

190 for i"l to n 

200 P-P«(l+r/b)'H> 

210 Print l;t«b(€);ii|ic;tiib(14);int(pitiie0+.5)/'100 
220 P"P+c 
230 next 1 



C»lcuUte» inter e»t *nd b«Unce for « »!(»tei»»tic savings Pro9r»m. 

Monthly dePotlt *? 10 
ComPoundlnS Periods Per year? 12 
Interest rite <xx.x>? 15 
No. of Ke«rs? 5 



End 




Current 


Ve«r 


Invested 


Balance 


1 


120 


139.29 


2 


240 


300.97 


3 


360 


488.65 


4 


480 


706.49 


5 


600 


959.35 
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CalcuUtM inttrtst tnd balance for « systematic savinst Pro9ram. 

Flonthlv deposit t? 10 
ConPoundlnS Periods Per year? 1 
Interest rate (xx.x>? 8.S 
No. of years? 8 



End 




Current 


Vear 


Invested 


Balance 


1 


120 


130.2 


2 


248 


271.47 


3 


360 


424.74 


4 


488 


591.04 


5 


680 


771.48 


6 


728 


967.26 


7 


848 


1179.68 


8 


960 


1418.15 



Calculates Interest and balance for a systematic savlnss Pro9ram. 

Monthly dePoslt *7 18 
Compound ln9 Periods Per year? 365 
Interest rate (xx.x>? 8.5 
No. of years? 8 



End 




Current 


Vear 


Invested 


Balance 


1 


120 


130.64 


2 


240 


272.88 


3 


360 


427.73 


4 


480 


596.32 


5 


600 


779.86 


6 


720 


979.68 


7 


848 


1197.23 


8 


960 


1434.88 



Try to combine the things that you have learned from the programs in 
these sections to write a completely generalized program for systematic sav- 
ings. It should accept the following information as input: 

• Initial deposit 

• Frequency of periodic deposits 

• Amount of each deposit 

• Interest rate 

• Interest compounding frequency 

• Length of time of savings program 

Your program should produce output in tabular form showing years (or 
other time periods), amount invested without interest, total amount with in- 
terest, and the interest alone. 
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Loan Payments 



Although saving money in a systematic way is a noble goal, many people 
frequently find themselves talking to a different bank officer, namely the one 
in charge of loans. 

This program will calculate the payments for a loan for a period of one 
year or longer. The program asks you to enter the key facets of the loan in 
question: amount borrowed, annual rate of interest, interval between pay- 
ments and term of the loan in years. 

The sample run shows a relatively short-term loan (2 years) of $3000 at a 
bargain basement interest rate of 8.5%. The monthly payment of the loan is 
computed to be $136.37 and the total interest $272.79. Why is the total in- 
terest not $255 (8.5% of $3000)? Instead it is 9.09% Is this a mistake in the 
program? 

Try the program for some longer term loans at real world interest rates. 
For example, run it for an automobile loan of $8000 at 12% over a 5-year 
period. Perhaps you can see from this that systematically saving your money 
and paying cash for an item makes more sense than time payments. 

Run the program to calculate the mortgage payments on a $150,000 house 
with a $40,000 down payment. Try it with a 16% interest rate stretched over 
a 30-year period. Ouch! Look at all that interest! 

At the heart of this program are Statements 180-210. Can you see what is 
happening in these calculations? 
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5 Print chr*<14>:rcn switch louitr Ctte 

10 Print "C«lcuUtt» PMnient schedule on « loan. " 

20 Print 

30 inPut "Amount borrowed";* 

40 inPut "Interest r«te <xx.x)";i 

30 inPut "Pwment interval (months) ";P 

60 InPut "Term <Ke«rs)";K 

70 input "Output table <!> or totals onli* <2>";v 
80 Print 

150 if v-2 then 180 

1£0 Print "Period Interest Principal" 

170 Print " Due Due " 

180 z-<s(#12>/'P 

190 k-<l«<P/'12»^100 

200 e»a#k,'<l-l/'<l+k>tz> 

210 e-int<e«100+. 5)^100 

220 c«a 

230 f-0 

240 dl-0 

230 tl-0 

260 tl-tl+1 

270 if tl>z th>n 380 

280 b-tl 

290 c"c-f 

300 d-ciKk 

310 f«e-d 

320 d>int(d«100-i-.5)/'100 

330 f-int<f#100+.3>/'100 

340 dl'dl+d 

330 if v-2 then 260 

360 Print b;tab<8>;d;tab<li);f 

370 9oto 260 

380 if v"l then 430 

390 dl"int<dl*100+.5)^100 

400 Print ''Principal *";a 

410 Print "Interest *";dl 

420 9oto 430 

430 Print " " 

440 Print "Total ";dl;tab<20);a:prlnt 

450 e3»int<<dl+a>*100+.!)/'100 

460 e6-e5/<(a#12VP> 

470 e6=inta00i|ie6+.5)/100 

430 Print "Total *";e5 

490 Print "Each Paanent t";e6 



Calculates Payment schedule on a loan. 

Amount borrowed? 3000 
Interest rate <xx.x)? 8.5 
Payment interval (months)? 1 
Tern (years)? 2 

Output table (1) or totals only (2)7 1 



Principal * 3800 



Interest « 272.79 
Total * 3272.79 
Each Payment « 136.37 
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Interest on Credit Purchases 



All too frequently, the rate of interest to be charged on a loan or credit pur- 
chase is hidden in the small type. After all, the bank or car dealer or finance 
company wants to convince you that you can afford that new car or home 
improvement or whatever. 

This program calculates the interest rate on a loan given the principal value 
of the loan, the number of payments, and the amount per payment. To make 
things even easier, the program will accept the cash purchase price of the arti- 
cle and the down payment and automatically compute the principal value of 
the loan. 

The sample run shows a rather unrealistic loan for an item costing $88.99. 
The down payment was $10.00 and the loan is over a period of 18 months; 
each monthly payment is $4.85. The program run shows that the actual in- 
terest rate is a modest 7.01%. 

Now try the program with a more realistic loan. A termite company in 
New Jersey advertises a total treatment for your house for only $200; just 
$29.95 down and 24 monthly payments of $1 1.95 per month. Is this the bar- 
gain that it seems? What is the annual rate of interest? 
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5 Print chr*< 14) : rem »wltch lowtr cam 

10 Print chr*(147)J"Int»rt»t r»tt of i Purchase on credit" 

20 inPut "Purchate Price";P 

30 input "Down P»Kiiitnt";el 

40 InPut "Number of Pt»mentt";n 

30 input "Payment* Per month" ;m 

60 input "Amount Per Payment"; a 

70 Pl-P-d 

80 t«a#n-Pl 

90 v-n/'<m«12> 

100 r«lnt<10000#t^Pl/'«+.5)/'100 

110 Print:print "Interest rate is";r;"X" 



Down Payment? Interest rate of a Purchase on credit 

Purchase Price? 88.99 

Soun Payment? 10 

Number of Payments? 18 

Payments Per month? 1 

Rmount Per Payment? 4.8S 

Interest rate is 7.01 y. 



Interest rate of a Purchase on credit 

Purchase Price? 1000 

Down Payment? 

Number of Payments? 24 

Payments Per month? 2 

Amount Per Payment? 50 

Interest rate is 20 '/. 
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Population Growth 



So far, all the compounding problems in this chapter have involved 
money — interest, savings, and loans. But many fascinating compounding 
problems do not involve money. Consider the following problem. 

In 1960, the population figures for the United States and Mexico were 180 
million and 85 million respectively. The annual population growth rate for 
the United States was 1.23% and for Mexico 2.23%. If these growth rates 
remain steady, in some distant year the population of Mexico will exceed that 
of the United States. In what year will this occur? 

The program uses the method of accumulating the principal amount rather 
than a formula, and applies it to both populations. The sample run reveals 
that the population of Mexico will overtake that of the U.S. in the not too 
distant future. 



5 Print chr*<14>:r«m switch lower ms* 

la u- 183009000 

20 m-8S000e00 

30 rl-1.0123 

40 r 2-' 1.0223 

50 »«19Se 

60 u"u«rl 

70 nt"nt*r2 

80 !*"K+1 

90 if i>i<u thtn 60 

100 Print "Vear 

110 Print t*b<9).: "Population" 

120 Print "U.S.R. ";u 

130 Print "Mexico ";iti 



Vear • 2037 

Population 
U.S.fl. 461406338 
Mexico 464464688 



But what if the Mexican people were dihgently trying to bring their 
increasing population under control and their annual growth rate was 
increasing by only 0.001% per year, compared to the increase of the U.S. 
growth rate of 0.01% per year. If that were the case, would the population of 
Mexico ever exceed that of the U.S.? 

Insert Statements 60, 70, and 85 into the first program, run it and find out 
the answer to the question. 

60 u»ij.*<1.01*rl) 

70 i»i=n*a.001«r2) 

85 if s>9999 then 110 
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Going back to the first program, run it with the population data from the 
U.S. (180 million) and California (15.7 million) from 1960. At that time, the 
annual growth rates were 1.23% and 3.7% respectively. Running the pro- 
gram will indicate the year that the population of California will exceed that 
of the United States. This, of course, is nonsense. Where is the discrepancy? 

This program might make more sense if it were restated to ask in what year 
the population of California would exceed that of the remainder of the U.S., if 
ever? 

Compounding can also be used to solve other kinds of population growth 
problems. Consider the bristleworm. The bristleworm can reproduce by split- 
ting itself into 24 segments, each of which grows a new head and tail. What is 
the maximum number of bristleworms that could be obtained in this fashion, 
starting with only one worm, after ten spHttings? Assuming a splitting occurs 
every 22 days, how many offspring will one worm produce in a year? When 
will the earth be overrun with bristleworms? 

Here is another problem for which the principle of compounding is useful. 
It takes nature about 500 years to produce one inch of topsoil. Many years 
ago, the United States h^d an average depth of almost nine inches of this 
good dirt, but as of 1975, the country was down to about six inches. This type 
of soil is necessary, of course, for growing food. 

Careless management of our soil causes about 1% per year to erode away; 
it is then lost forever. Once soil depth reaches three inches or less, it is impos- 
sible to grow crops on a large scale. Write a program to calculate the year in 
which the U.S. will have less the 3" of topsoil, assuming that it continues to 
erode away at 1 % per year. 

Will you be alive then? Will your children be aHve? Will anyone be alive? 
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Probability 



Statistics and probability are subjects that bring up all kinds of images. 
Some people who have not had a pleasant time in a required college statistics 
course keep as far away from the subject as possible. To other people, statis- 
tics is something with which devious manufacturers can mislead you about 
their products. 

Not long ago, one foreign automobile manufacturer boasted that 90% of 
their cars sold in the United States in the last seven years were still on the 
road. This sounds like excellent reliability. But consider the fact that this 
manufacturer was rapidly expanding in the U.S. market and 65% of the cars 
they had sold in the U.S. had been sold in just the previous two years. You 
would expect that all of these would still be running. 

If the manufacturer had sold 10% of their 7-year volume in Years 1 to 3, 
and most of these cars were not running, then the advertising claim loses 
much of its meaning. 

When approached in a logical, step-by-step manner, statistics and probabil- 
ity are not difficult. In fact, they can be a great deal of fun. 

Some of the programs use a formula while others simulate the event in 
which we are interested. The results are the same, but you may find that the 
simulations help you to understand exactly what is happening. 
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Pascal's Triangle — Calculated 



Pascal's Triangle is quite fascinating. As you can see from the portion of 
one reproduced below, each row is symmetrical. Each row also contains the 
coefficients for a binomial expansion. The sums across the ascending diag- 
onals form the Fibonacci sequence. The sums across the rows are all powers 
of two. Each row corresponds to the digits of a power of 1 1. Every element is 
the sum of the two above it. And, in case you care, all the elements in it are 
identities in combinatorial theory. 

The ways this marvelous triangle can be generated are as varied and in- 
teresting as its properties, though perhaps more difficult to figure out. Here is 
one way to use a computer to generate the triangle. 

Any element can be found be adding together the two elements immedi- 
ately above it. The program uses this principle to produce a triangle eight lev- 
els deep. Lines 10-30 set the rightmost diagonal to 1. Each element is stored 
in the two-dimensional variable P(R,C) with R denoting the row and C 
denoting the "column." The variable T (Line 50) simply leaves some blank 
spaces so the output resembles a triangle. The crux of the calculation is in 
Line 70 in which the value of each new element is calculated. 

There is another interesting way to calculate Pascal's Triangle in even 
fewer statements than the program here. It generates the triangle one element 
at a time and does not use any arrays or two-dimensional variables. Can you 
figure out how to write such a program? 



10 for c"l to 8 1 
20 i><e.e>«l 

38 ntxt c 1 1 

40 for r«l to 8 

90 t-20-r«2 1 2 1 

£0 for C"2 to r+1 

70 P<r+l,c>«P<riC)4*(r,c-l) 13 3 1 

80 If l>(r<c>-0 thtn 110 

90 Print ttb<t>;P(:r,c>; 1 4 S 4 1 

100 t»t+4 

110 ntxt c 1 S 10 10 5 1 

120 Print 

130 Print 1 6 15 20 15 6 1 

140 ntxt r 
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Pascal's Triangle — by Probability 



This program simulates the dropping of balls through a triangular array 
shown below. At each level, a ball is equally likely to fall either to the left or 
right. At the bottom of the array a cup is placed at each end point; these cups 
collect the bails. After each group of balls has been dropped, the number of 
balls in each cup is tallied and displayed. 



/\ 
/\/\ 
/\/\^ 
/\/\/\/\ 



u u u u u u 



S Print chr«(i4>:rtm tultch loutr cut 
le Print "Puctt't trtlnStt by ProbablUtv" 
20 inPut "BtUt to droP";m 
30 inPut "No. of l«ve1»";k 
40 Print "Petition B*Ut" 
50 for n"l to m 
£0 t-0 

70 for 1"1 to k 
80 if rnd(l».5 thtn 100 
90 t-t+1 
100 ntyt 1 
110 b<t+l)-b(t+l>+l 
120 n«xt n 
130 for l-l to k+1 
140 Print 1it«b(12>;b<1> 
130 b(l)-0 
160 ntxt 1 



P««c*r» tminSl* bu Probability 
B«Ut to drop 1000 
No. of Uvttt 2 
Petition B«Ut 

1 249 

2 St£ 

3 233 



Pascal 't tralnait bv Probability 
Ballt to drop 1000 
No. of 1*vt1t 4 
Petition Ballt 



£4 

241 

390 

241 

£4 
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Pji»c*1'» tr*ln9l» by ProbabllitK Pascal's tralnsle ba ProbabilltK 

Balls to drop 160 Balls to droP 640 

No. of Itvtls 4 No. of l»w«ls 6 

Position Balls Position Balls 

19 19 

2 45 2 70 

3 56 3 142 

4 40 4 208 

5 10 ! 134 

6 51 

7 6 



The number of balls landing in the various cups at each level, when re- 
duced to the lowest common denominator, should approximate the numbers 
in Pascal's Triangle. Do you know why they should? 

Try a few runs of the program with a different number of balls. Do the 
numbers obtained really approximate those in Pascal's Triangle? If, at each 
dividing point in the array, every other ball went in the opposite direction, 
then Pascal's Triangle would be produced. Try running the program with 4 
balls at Level 2; with 8 balls at Level 3; and with 16 balls at Level 4. Do your 
results look like those of the previous program which calculated the triangle? 
Does the approximation come closer as the number of balls is increased? 

How can you determine how close your results from this program are to 
the exact value of Pascal's Triangle? One way is to take the number of balls 
that dropped into each cup on a given level and divide that by the total num- 
ber of balls divided by the theoretical sum of the row. So, in the sample run at 
Level 4, you would divide the balls in each cup by 62.5 (1000/16). You then 
compare that to the "correct" number and compute the percent difference. 

Here is the result of this procedure for Level 4 (actually the fifth row) in 
the sample run: 



Cup Balls —62.5 Correct Deviation 

1 64 1.02 1 2.00% 

2 241 3.86 4 3.50% 

3 390 6.24 6 4.00% 

4 241 3.86 4 3.50% 

5 64 1.-2 1 2.00% 
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Common Birthdays 



Here is an interesting little problem in probability. In a group of ten people 
selected at random, what is the probability that any of them will share the 
same birthday? How about a group of 20 people? Of 50 people? 

Conversely, how many people would you need in a group such that there is 
a 50% probabiUty that at least two of them have the same birthday? How 
many people would be needed for a 90% probability of an overlap? 

Try to answer these questions, either by guess or by calculation, before you 
look at the output from the program. 

This program provides a painless introduction to the world of statistics. 
The calculation is actually quite trivial. The probabiUty that any person in a 
group has a birthday on January 1 is 1/365. If our group has only two people 
in it, the probability that both of them have a birthday on January 1 is 1/365 
times 1/365, or a very small number indeed. The probability that they have a 
common birthday is 365 times the very small number just obtained, or about 
0.27%. 

However, if there are three people in the group, the probability goes up 
sUghtly. Call the three people Betsy, Ken, and Larry. From the reasoning 
above, we know that there is a 0.27% probability that Betsy and Ken have 
the same birthday; also a 0.27 % probability that Betsy and Larry share a 
birthday; and finally, a 0.27% chance that Larry and Ken have a common 
birthday. Hence, the total probability for a group of three people is three 
times the probability for just two people. 

A group of four increases the probability over that of two people by a fac- 
tor of six, five people by a factor of 10, six people by 15, seven people by 21, 
and so on. There is a progression here, but the probability can also be cal- 
culated by the formula: 



Can you see why this formula produces the same result as the description 
in words and associated progression above? 

Consider all the presidents of the United States. (How many have there 
been to date?) Two of them, James Polk and Warren Harding, were born on 
November 2. Is this to be expected given the size of the group? 

Since birthdays can be predicted, at least statistically, it ought to be pos- 
sible to predict deaths as well. It is interesting too that John Adams, James 
Monroe, and Thomas Jefferson all died on July 4th. Millard Fillmore and 
William Taft both died on March 8. What is the probability of that set of 
events? 
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5 Print chr*<14>:ren switch tower c«*t 

10 Print "People Sane birthdas ProbabiUtK" 

20 <i-364/'3S5 

30 for n"2 to 40 

40 P-ie0W(l-<1> 

50 Print n;tab(20);int(P«100+.5>/'100;"X" 

60 i-qmoss-nJ/ssi! 

70 next n 



PeoPle Same blrthdaK ProbabiUt! 

2 .27 Ji 

3 .82 ^ 

4 1.64 X 

5 2.71 K 

6 4.05 K 

7 5.62 X 

8 7.43 X 
5 9.46 X 

10 11.69 X 

11 14.11 X 

12 16.7 X 

13 19.44 X 

14 22.31 X 

15 25.29 X 

16 28.36 X 

17 31.5 X 

18 34.69 X 

19 37.91 X 

20 41.14 X 





44.37 X 


22 


47.57 X 


23 


50.73 X 


24 


53.83 X 


25 


56.97 X 


26 


59.82 X 


27 


62.69 X 


28' 


65.45 X 


29 


68.1 X 


30 


70.63 X 


31 


73.05 X 


32 


75.33 X 


33 


77.5 X 


34 


79.53 X 


35 


81.44 X 


36 


83.22 X 


37 


84.87 X 


38 


86.41 X 


39 


87.82 X 


40 


89.12 X 



Now it is your turn to write a program. Here is a game to be played among 
your friends. You make a bet that they have similar preferences in colors. 
Each person in the group puts up one penny and you, because you are so sure 
of yourself (and probability), put up an amount in cents equal to the number 
of people in the group. If you win, you get to keep all of the money; if your 
friends win, each of them gets double his original bet, or two cents. 

Then each person selects a color (from a larger list than there are people, of 
course). If any two have picked the same color, you win; if all have picked 
different colors, they win. 

Since you would like to win, you need to know how many colors should be 
on the list for different size groups to give you a better than 50-50 chance of 
winning. You might like the probability to be up around 70% or so. You can 
produce the necessary table with a five-line program. Go to it! 
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Coins in a Poclcet 



The next two programs were originally written by Glenda Lappan and 
M.J. Winter and appeared in Creative Computing magazine. They are marvel- 
ous simulations for illustrating various aspects of probability. 

Coins in a Pocket is a simulation of the following situation. A newspaper 
costs 5 cents. A customer has 5 pennies and a dime in his pocket and offers to 
pay for his paper by letting you, the vendor, select at random two of the six 
coins. If you and this customer repeat this procedure for the 20 working days 
of a month, how much more or less than $1.00 (20 days x 5 cents) are you 
likely to have collected? 

The program below solves this little problem. The first random coin is se- 
lected from a group of six (Line 50 in which X = INT(RND*6)). The value 
of X can be 0, 1, 2, 3, 4, or 5. If it is 0, we assume the dime was selected and a 
second pick is not made, because it will surely be a penny. Thus 1 1 cents is 
added to the running total in Line 1 10. 

If the first coin is a penny (X = 1, 2, 3, 4, or 5) then we make a second pick 
from the five remaining coins. Again, if the dime is chosen (Y = 0), 1 1 cents 
is added to the total; otherwise 2 cents is added. 

As you can see from the sample runs, after a great number of trials, the 
average amount of money collected each day seems to be close to 5 cents. If 
you would like to convince yourself that the answer is, in fact, 5 cents, con- 
sider the following. Assign a letter to each coin, A - F. Let A be the dime and 
B through F be the pennies. Since all combinations are equally probable, here 
are all the possible combinations: 

AB 

AC EC 

AD BD CD 

AE BE CE DE 

AF BF CF DF EF 
There are five combinations with a dime (5x11 cents) and ten combina- 
tions of only pennies (10 x 2 cents). Add it up and you have 55 cents plus 20 
cents divided by a total of 15 combinations which equals an average value of 
75/15 or 5 cents. 



88 



5 Print chr»<14)irt» loutr cii«t 

10 Print "SimuUttt ttkin9 2 coins nt r«ncion from « dint tnd 5 Ptnnin" 
20 inPut "No, trUU"Jn 
30 u"0 

40 for k'l to n 
S0 x-lnt(rnd<l>i«6) 
£0 If x-0 thtn 110 
70 »-int<rnd<l>«6> 
S0 if y>0 thtn 110 
90 u"u+2 
100 goto 120 
110 u-u+11 
120 ntxt k 

130 Print "flvtr«9t vtlut it";u/'n;"ctnt»." 



Simulatts t«klne 2 coins «t random from a dlmt and 5 Ptnnits 

No. trials? 100 

RvtraSt valut is 4,16 ctnts. 



Simutatts taklnS 2 coins at random from a dlmt and 5 Ptnnits 

No. trials? 1000 

flvtraflt valut is 4.7 ctnts. 
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Baseball Cards 



In statistics and probability, it is frequently necessary to predict the num- 
ber of trials on average until one is successful. For example, if you are trying 
to roll a four with one die, on average how many tries will it take until you do 
so? 

Since a die has six sides, the probability of rolling any number between 1 
and 6 is p = Vg. Thus, on any given roll, you have a Yg chance of rolling a 
four and chance of rolling something else. This failure is designated q. So 
to be successful in rolling a four, we have a Ve chance on the first roll. On two 
rolls, the probability is two trials x the probability of one failure, then success 
(2 X Ye xVe). On the third roll, the probability of success is three trials x the 
probability of two failures followed by one success (3 x % x % x Ve). 

Continuing this reasoning leads to the formula for the expected number of 
trials to success: 

E=lp + 2qxp + 3qxp + 4qxp + ... 

This series can be reduced and solved for E which leads to: 



l-q p 

So now the answer to the original problem can be calculated. The expected 
number of trials to roll a four is 1/p = 6 tries until success. 

But there is another way to approach this type of problem with the assis- 
tance of the computer. Consider the problem. Assume there are ten different 
prizes in Crunchies cereal boxes. How many boxes of cereal would you have 
to buy to obtain the complete set? The formula above can be expanded to 
solve this problem: 

For a value of ten, you can solve this formula by hand. However, let's say 
you would like to solve this problem for baseball cards as found in packs of 
bubble gum. If there are 50 cards in a complete set, how many packs of gum 
must you buy, on average, to get a complete set. Write a computer program 
using the general formula above to solve for a set of any size. The answer for 
50 cards is 224.96 packs of gum; how many would you have to buy for a set of 
100 cards? 
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There is another way of approaching this problem. This method is similar 
to that used to randomly select coins out of the pocket. In this case, the ran- 
dom group is set equal to the total number of cards; this value is accepted as 
input in Statement 20. Each purchase is set equal to a random number be- 
tween 1 and N (Statement 120). This card is then put into its proper place in 
the collection (Statement 130). However, if there is already a card there from 
a previous purchase, we simply increment the purchase counter (Statement 
160) but do not get any closer to obtaining a complete set. After each pur- 
chase, we test to see if the set is complete (Statement 170), otherwise we go on 
buying more packs of bubble gum. 

When the set is complete, the number of packs of gum are tallied up and 
printed out. After a set of trials, the average is computed. This averaging 
value should be reasonably close to the value obtained by the formula al- 
though it will take a great number of trials before the two numbers are within 
1 % of each other. 



5 Print chr*(i4>:rtm switch loucr cmi 

10 Print "SlmuUttt buvlnQ Sun ulth bMtbtU ctrdt." 

20 InPut "Ccrds In Mrltt";n 

30 inPut "No. ilmlUr c«rd*";k 

40 dim c(100) 

50 Print "TrUl P*ck»" 
60 «-0 

70 for 1-1 to k 
90 da0 

90 for J"l to n 
100 c<J>-0 
110 ntxt i 

120 x-int(rnd<l>«n>')-i 
130 c(x>ac(xHl 
140 if e<x>»l thtn 1S0 
150 Soto 120 
160 d-d^l 

170 if d-n tKwn 190 
180 eoto 120 
190 t-0 

200 for J'l to n 
210 t-t+c<J> 
220 ntxt J 

230 Print i;t«b(ll);t 
240 ft*t 
250 ntxt i 

2fi0 Print "nvtr««t "it/k 
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SlmuUtM buvlnS 9um with b«stb«ll 
card*. 

Cards lit ttritt? 10 
No. simiUr cards? S 
Trial Packs 

1 36 

2 2S 

3 21 

4 20 
3 19 

flvtraSt 24.2 



If you run this second program for a set of 100 cards, it will sometimes run 
for a very long time before arriving at the answer. Try it with a set of 500 
cards as are used by some real bubble gum manufacturers. You can go have 
your dinner while the program computes just the first value. Be sure to 
change the dimension statement in Line 40 to C(500). 
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System Reliability 



As more and more people in the world come to depend upon mechanical 
and electronic devices in a myriad of different ways, it is important that these 
technological devices continue to function. 

Some years ago, the military came up with a measure that could be applied 
to all kinds of systems, big and small, to measure reliability. It is called 
"mean time between failures." What this means is the length of time, on av- 
erage, between breakdowns. For a tank, this may be 100 hours, while for a 
spacecraft, the MTBF must be considerably longer than the planned mission. 

As we saw in an earlier section, it is frequently desirable to break down a 
large problem into smaller subproblems. To calculate the MTBF for a space- 
craft would be quite impossible. Instead, it is necessary to start with smaller 
systems and build up to the whole. 

Consider one of the electrical subsystems of a spacecraft. It uses five 
components as shown below, two parallel systems A and B, arranged in se- 
ries. The parallel subsystems are said to be "redundant." This is one method 
of increasing reliability since the system will continue to work if at least one 
of the components works. 




If the manufacturer of the components stated that each one has a 60% 
probability of lasting 1000 hours, what is the chance that the entire system 
will last 1000 hours? Take a guess before reading on to the solution below. Is 
your guess greater than 60% or less? Why? 

The program below is a simulation of this system. Remember in subsystem 
A, it will continue to work if any of the three components works and in B if 
either of the two components works. However, the system will fail if all three 
of the A components or both of the B components fail. 

In Statement 30, the program is set up to make 500 trials of the system. 
Statement 50 selects a random integer between 1 and 10 for each component 
for each trial. If the value of this integer is 6 or under, the component is okay 
(60% probabiUty of working at the end of 100 hours). If it is 7, 8, 9, or 10, 
this means the component has failed, and the program prints "Bad!" 
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5 Print chrt<l4):rem switch lower c»s» 

10 Print "SinruUtM u»» of « circuit for 1000 hr«":prlnt 

20 t'0 

30 for i"l to 300 

40 for J"l to 5 

30 c(J>»lnt<rncl<l)«10)+l 

£0 If c(J)>S then 90 

70 c<J>-l 

80 9oto 100 

90 c<J>«0 

100 next i 

110 if c<l>"l or c<2)«l or c<3)"l then 140 
120 Print "Btd!"; 
130 9oto 180 

140 if c<4)«l or c<5)»l then 170 
150 Print "Bad!"; 
160 goto 180 
170 »"s+l 
180 next 1 

190 Prlnt:print "RelUblUta is aPProx. lnt<10000#<*/'500)+.5V100; "K" 



Simulate* use of a circuit for 1008 hrs 

Bad ! Bad ! Bad ! Bad ! Bad I Bad ! Bad I Bad ! Bad ! Bad ! 
Bad ! Bad ! Bad ! B«d ! Bad ! Bad ! Bad ! Bad j Bad ! Bad ! 
Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! 
Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! 
Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! 
Bad ! Bad ! Bad ! Bad ! Bad ! Bad I Bad ! Bad ! Bad ! Bad ! 
Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Sad ! Bad ! Bad ! Bad ! 
Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! 
Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! Bad ! 
Bad! 

ReliabiUta is aPPrex. 81.8 V. 



Statement 1 10 tests whether subsystem A works and Statement 140 checks 
out subsystem B. If both work, a success is counted in Statement 170. 

There are, of course, shorter ways to do this. For example, all the individ- 
ual tests (Statements 110 to 140) can be replaced by one overall test (new 
Statement 110). 



110 If <c<l)+cC2)+c<3))«<c<4>+c(3)>-0 then 140 
130 s»s+l 
140 next I 

RellabUits Is aPProx. 78.9 •/. 



Now it is your turn to use what you have learned in this chapter and write 
a program to determine the mean time between failures of the system above. 



94 



7 

Geometry and Calculus 



In this chapter several of the problem solving approaches from earlier 
chapters are used to solve geometric problems. You will find that many prob- 
lems can be solved in a variety of different ways — ^by applying a formula, by 
trial and error, by successive approximations, and, in some cases, by common 
sense. Perhaps the most important thing to learn from these problems and 
programs is how to analyze a problem to reach the solution quickly and 
painlessly. 
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Crossed and Slipping Ladders 



Here is a simple problem of a slipping ladder. A ladder 25 feet long is 
placed so its foot is 7 feet from the base of a building. The base of the ladder 
slipped on some loose gravel so that the top is 4 feet lower than where it was 
to start. How far did the foot of the ladder slip? 




7 Y 



The diagrams show the two positions of the ladder. By the Pythagorean 
theorem we know that a^ + b^ = c^, hence, the equations needed to solve the 

problem are: 

X = \/252 - V 



y = VIS^ - (x-4)= 

and the amount of shppage of the base is z = y - 7. These three equations are 
put into a computer program which quickly calculates an answer of 8 feet. 



10 Print "SllPPlnfl L«del»r" 
20 C'2S 
30 b-7 

40 X"»0r<ct2-b12> 
50 i(«»<)r<ct2-<x-4)t2> 
60 z"K-b 

70 Print "Bute »UPP«el";z;"ft" 



fllPPlnO Udder 
»»e slipped 7.99999999 ft 



While the aritmetic in the above problem is not particularly messy, it is still 
no great joy to solve by hand. However, the computer is just as happy to do 
the problem with really messy dimensions, say a ladder length of 27.83 feet 
and a distance from the wall of 7.62 feet. 

20 c»27.83 SllPPlnfl Ladder 

30 b-7. €2 Eiise tUPPtd 8.38612 ft 
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Let's consider an old problem found in many classic collections. Two lad- 
ders, one 25 feet long and the other 35 feet long lean against buildings on op- 
posite sides of an alley as shown below. The point at which the ladders cross 
is 12 feet above the ground. How wide is the alley? 




By using similar triangles twice we find that 

12/a + 12/b = 1 
Then, by applying the Pythagorean theorem and reducing, we obtain: 



a^ -b^ 



600 



Using one of the methods described in the Problem Solving chapter, you 
can solve these two simultaneous equations. Or they can be combined into 
one equation in which z, the width o f the alley is 

z = -/35^ - a' 

By eUminating b, the following fouth degree equation is obtained 
a* - 24a^ - 600a^ + 14400a - 86400 = 

Again, this can be solved using one of the methods in the chapter on Prob- 
lem Solving. This is the traditional approach, but there is another. The solu- 
tion is the intersection of the curves described by the two original equations. 
By using successive approximations, say in steps of 1 for a, you could solve 
for bj and \ in the following restated original equations 



12a 



12 



b^ = a/600 - a^ 

When bj falls below b^, reduce the step to 0. 1 to obtain a closer approxima- 
tion. By continuing this method of successive approximations, it is possible to 
obtain a very accurate solution. 
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Is there another approach? Yes, there is and it also avoids the quartic equa- 
tion. It uses the original equations in a trial and error procedure as described 
in the chapter on Sets and Repetitive Trials. See if you can write a program 
using this approach. 

Here is another classic problem for you to solve in any way that you like. 
What is the longest ladder that can be carried in a horizontal position around 
the corner made where a 12-foot wide alley meets one that is 8 feet wide. The 
diagram shows the problem. 
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Distance Between Coordinate Points 



This program solves for the distance between any two points in three- 
dimensional space defined by their x, y, and z coordinates. 

The formula for the distance between two points in three-dimensional 
space is: 

d = (X, - x,)^ + (y, - y,)=^ + (z, - Zj)=' 

It would be desirable for the program to be able to solve any problem of 
this kind. One approach would be to use input statements to accept the point 
coordinates. Another is to use a data statement to define the points. Then 
only this one statement has to be changed for a new problem or set of 
problems. 

The program is written to calculate the distance between three sets of 
points. The points used were: 

0, 0, and 3, 4, 5 

3.5, -4.7, 6.2 and -0.9, 3.0, 4.4 

67, 36, 82 and 54, 25, 90 

This calculation is used extensively in aerospace navigation. Can you deter- 
mine the angle or "compass heading" of the resulting flight path? Better start 
with just two dimensions. 

5 Print chr*<14) -rtn twitch lowtr c*«t 
10 Print chr*<147) 

20 Print "C«lcuUtcs dlttunc* bctuttn Points In 3-II sPact" 
30 Print 

40 Print "Coord In* te» Dl»t«nc»" 
50 rtJid i,b,c,d,eif 
60 l-s<lr<<«-d>t2+(b-»>12+<c-f>1«> 
70 Print *;b;c 

80 print d;t;f;t«b<18);lnt<l«1000+.5>/'1000 

90 Print 

100 9oto 40 

lie d«ta 0.0.0.3/4.5 

120 d«t« 3. 5. -4. 7. 6. 2. -.9. 3, 4. 4 

130 d«t* 67.36.82.54.25.90 

C*lcuUt«s distance between Points In 3-D sPace 
Coordinates Distance 



3 4 5 7.071 

Coordinates Distance 

3.5 -4.7 6.2 
-.9 3 4.4 9.049 

Coordinates Distance 
67 36 82 

54 25 90 18.815 
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Area — by Calculation 



It is a simple matter to calculate the area of regular geometric figures by 
using the usual formulae. However, the problem becomes more difficult when 
it is necessary to calculate the area of two combined regular shapes or the 
area of irregular shapes. This program shows the method of analysis for solv- 
ing a problem of the first type, while the next program demonstrates four 
methods of dealing with irregular areas. 

The problem is to solve for the shaded area of the figure below for any 
value of the radius, R. 




The first step is to recall the formulae for calculating the area of a circle 
and square: 

A (circle) = pi * 
A (square) = Side^ or (2 * R)^ 

The difference in area between a square and a circle inscribed within its 
borders is: 

A (difference) = A (square) - A (circle) 
and the area of one corner is the difference divided by 4. The program below 
will calculate the area for any radius. 

Now it is your turn. Write a program to calculate the difference in area 
between a circle and square in which the square is inscribed within the circle. 
Now, change your program to calculate the difference in area for a triangle, a 
hexagon and an octagon inscribed within a circle. 
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5 Print chr*<14> :rtn »witch lower c«sc 

le Print "flrM bttkiMn circle inscribed in square" 

SB input "R«eiiu«";r 

30 c-3. 14139«rt2 

40 >-(2ilir)'r2 

SB d»»-c 

68 HliiM 

70 Print "Tr»PP»d «re»";d 
80 Print "One corner" ;dl 



flret between circle intcribed in Muare 
Rtdiut 10 

TrtPPed irea 83.8410001 
One corner 21.46023 

nre« between circle inscribed in square 
Radius 1 

Trapped area .838409999 
One corner .2146023 



Extend your program to calculate the difference in area for any regular fig- 
ure inscribed within a circle of radius 1.0. Set up a table or results as follows: 
Number of sides Difference in Area 

3 
4 
5 
6 
7 
8 

What can you conclude from these results? Do these areas follow in some 
sort of progression? 
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Area^by Integration 



In many cases it is necessary to determine the area of an irregular figure or 
the area under a curve where an exact formula is not available. The approach 
most commonly used is to divide up the enclosed area into small regu.i'-ly 
shaped pieces and sum up the areas of all of these pieces. 

The easiest shape to use in these calculations is a rectangle. A group of 
rectangles can either be inscribed within the irregular figure or curve, or 
circumscribed around it. The first two diagrams show these two methods be- 
ing used to find the area under a curve. A third method is to use trapezoids 
which, depending upon the direction of curvature, will either be inscribed or 
circumscribed automatically. 




Inscribed Circumscribed Trapezoid 

Rectangle Rectangle 
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A fourth method, known as Simpson's Rule; essentially fits a series of 
parabolas between the points of the curve and calculates the average area. It 
requires that the area be divided into an even number of parallel slits. Let us 
call the total number of divisions 2m which are h distance apart. The point on 
the curve where the first line intersects is yo, the second yi, and so on until 
y2in- The area is then given by the formula: 

A= y3h[(yo + y2m) + 4(y, + yj + ... 

+ y2n,-l) + 2(y2+ y4 + ... y2m-2)] 

The area as calculated by Simpson's Rule converges extremely quickly 
compared to the other methods. Nevertheless, as the number of intervals in- 
creases, all the methods approach the same limit. As might be expected, the 
trapezoidal approach converges more quickly than either of the approaches 
using rectangles. Compute the average of the two methods using rectangles. 
What do you get? Does this suggest another method? 

The methods used in this program involve the calculus. And you thought 
the calculus was difficult! Now you know otherwise. 

5 Print chr*<14) :rfin twitch loutr east 
10 print chr«<l<»7> 

20 Print "CoinPutt* trt undtr « curvt bn 4 ntthod* of int»9r«tion." 
40 Print 

50 Print "Enttr »t«rt «nd end v*luM for x (*n*ll ont first)" 
£0 inPut ii<b 

70 Print "R«n9t of x-";«;"/';b 

80 Print: Print "No. intervals 

90 Print " Inscribed Rtcts TraPtioids" 

100: Print "Circumscribed Rects P«r*boUs" 

110; n— 2 

120 s-0 

130 def fn«<x>x19 
140 n"iii+3 

130 for n"m to n+2 
160 c*s 
170 ^-0 
180 P-0 
190 d«2tn 

209 b'<b-t.)/6 

210 for 1-0 to d-1 
220 x-A+iDih 

230 P"P+h«fn»<x> 

240 q-q+hHifnaCx+h) 

250 next i 

260 t-<P+1>^2 

270 u«fn«<«)+fn»<x+h) 

280 for >2 to <d-2) steP 2 

290 a"u'+2«fn«<«+J«h> 

300 neyt j • 

310 V-0' 

320 for k«l to, <d-l) step 2 
330 viv+4#fn*<«+k#h> 
340 next k 
330 S"<u+v)«<h/'3> 
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360 Print:Print d:Prlnt:Prlnt P;t*ba2>;t:Print q;tab<12).;s 

370 next Ti 

380 If d<64 then 140 

390 if abs«c-s>/'<<c+«)/2>».0001 then 140 
400 stop 



Computes «ri under « curve by 4 nethodt of lnte9r«tion. 

Enter st«rt «nd end values for x (sm«U one first) 

? 1 
?? 10 

Ranee of x> 1 , 10 
No. Intervals 

Inscribed Rects Trapezoids 
Circumscribed Rrcts Parabolas 



2 

753.187501 
5248.6875 



30d0.9375 
5499.75001 



4 

1501.17188 
3748.92188 



2625.04688 
2499.75 



8 

1969. 13672 
3093.01172 



2531.07422 
2499.75 



16 

2226.61231 
2788.54981 



2507.58106 
2499.75 



32 

2361.22339 
2642. 19214 



2501.70776 
2499.75 



64 

2429.99726 
2570.48163 



2500.23944 
2499.75 
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8 

Science 



Using techniques and approaches presented in the previous chapters, this 
chapter contains five programs in the area of science. One uses a formula to 
solve simple gas problems, two are drill exercises on the gas laws of Boyle and 
Charles, and the last two are simulations. You will see that the simulations 
draw upon many previous techniques such as progressions and repetitive 
calculations. 
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Gas Volume 



Here is a simple program to produce a not-so-simple table of values for gas 
volumes. 

The volume of a gas varies directly with the absolute temperature T (Kel- 
vin) and inversely with the pressure P. If a certain quantity of gas occupies 
500 cubic feet at a pressure of 53 pounds per square foot and an absolute tem- 
perature of 500 degrees, what volume will it occupy at 600 degrees absolute 
temperature and pressures from 100 to 1000 pounds per square foot in in- 
crements of 50 pounds? 

The original conditions are used to solve for the constant K in Statement 
50 (K = V*P/T). Then new volumes are computed for varying pressures 
with T equal to 600 degrees. The formula used is V = K*T/P. 

In the second part of the program, Lines 70 to 100 are replaced to produce 
a plot of the gas volume for the various pressures. 

How would you modify the program to deal with a more general case (i.e., 
other gasses and different temperatures)? Second, can you write a program 
that produces a table of values for a gas at different pressures and 
temperatures? 



3 Print chr*<14):r»n twitch lowtr ctst 

10 Print "Om volunt* «t different Prttiurti" : Print 

20 v"500 

30 P-53 

40 t-S00 

50 kmv$P/t 

60 inPut "TtioPtrtturt (l<)";ti 
70 Print "PrMture Volun*" 

80 for P«100 to 1000 »ttP 50 OM votumts «t different Pressures 

90 v"k*tl/'P 

100 Print P;t*b(10>;v Temperature (k>? £00 



Pressure 


Volume 


100 


318 


150 


212 


200 


159 


250 


127.2 


300 


10£ 


350 


90.8371429 


400 


79.5 


450 


70. £££££€7 


500 


S3.£ 


550 


57.8181818 


£00 


53 


£50 


48. 92307fi9 


700 


45.4285715 


750 


42.4 


800 
850 


39.75 


37.4117fi47 


m 


35.3333333 


950 


33.4736842 


1000 


31.8 
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5 Print chr*<14) :r«iit switch lower cm* 

10 Print "0«» voluoiM *t dlfftrent Prt»»ur««" : Print 

20 v-see 

30 P-S3 
40 t-S00 
30 k>v«P/'t 

60 input "TtnP«r»turt <k)'";tl 
68 Print 

70 Print "Pr«*»urt Volume Plot" 
80 for P-100 to 1000 step 50 
90 w-k#tl/'P 

100 Print P;t«b(4+v/'17);"«" 
110 next P 



Gas voiuMcs *t different prassures 



TMHierature <k>? G8fl 




95» » 
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Charles' Law Drill 



In the previous program, gas volumes were calculated using Charles' Law 
and Boyle's Law. Here is a drill and practice program (remember Chapter 1?) 
that produces problems relating the volume and temperature of a gas. When 
pressure is constant, the volume/temperature relationships can be stated as 
follows: 

Vq ^ V, 
To T, 

The program presents four problems, one to solve for each of the four vari- 
ables in the equation above. 

How can you make the program more efficient? More interesting? 



S Print chr»<14) :rtiii twitch lou«r c«*t 

10 Print "Ch«rlM' Liw drlU" 

20 Print "Velum* In milliliter*" 

30 Print "TtmP*r«tur» In d»9r«»* Ktlvln" 

40 Print 

50 v-lnt<rnd<l>»50+50>il00 
S0 ul-int(rncKl>«S0+S0>«l00 
70 t-lnt<rncl<l>«123+l25) 
80 tl-ltit<rncl(l)*123+123) 
90 c-c+1 

100 on c «oto ll0.1S0;190,230>27e 
110 v-0 

120 Print: Print "Solve for w uhtn" 
130 <ll»vl«t/'tl 
140 Ooto 280 
1S0 vl-O 

160 Print " What It vl Slvtn" 
170 <ll-v*tl^t 
180 goto 280 
190 t-0 

200 Print: Print "CtlcuUtt v*lut of T" 
210 <ll«v*tl/'vl 
220 Ooto 280 
230 ti>e 

240 Print: Print "Solve for tl Slven" 
230 <ll-t#vlXv 
260 ooto 2S0 
270 ttoP 

280 Print "v- "; v;t«b<ll); "t-";t 
290 Print "vl" ";vi;t*b<n); "ti-";ti 
300 InPut "Vour «ntuer";q 
310 Print "Correct vilue";fll 
320 Print 
330 9oto50 
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Ch»rl»»' L«w drill 
Volume in milllllttrs 
TtmPtrnturt In d»9i"et» Ktlvln 



Solve for v when 
V- t- 217 

vl> £000 tl- 196 

Vour jintuer? 114418 
Correct vuliie 6642.95715 

Mhat i« vl given 
V" S000 t" 153 

vl- tl- 243 

Vour an«uer7 523821 
Correct vulue 7941.17647 



CclcuUte viluc of T 
V- 7100 t" B 

vl- 9800 tl- 162 

Vour tnswer? 18494 
Correct vjilue 117.367347 



Solve for tl Slven 
V- 9500 t- 222 

vl- 8300 tl- 

Vour «n»uer? 75000 
Correct v»lue 193.957895 
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Boyle's Law Drill 



Boyle's Law describes the behavior of gases under ideal conditions when 
pressure and volume are varied. When the temperature is constant, Boyle 
found that: 

Po * Vo = P, * V, 
Pressure is normally measured in centimeters of mercury while volume is 
measured in millilitres. As with the drill on Charles' Law, this program 
presents four problems, one to solve for each of the four variables in the equa- 
tion above. 

Unlike many other drill and practice programs, these two do not compare 
your answer with the correct one. Instead, they leave that up to you to do. Is 
this desirable? Why or why not? If you feel it is undesirable, change the pro- 
grams so they do compare answers and calculate a score. 



3 Print chr»(14>:rei>i switch louitr c«tt 
10 Print "loKlt't Uw drill" 
20 Print "Volumt In MlUlUttrt" 
30 Print "PrMsure In cmmer"; 
40 Print 

50 w-lnt<rnd<O#50+5a)«100 
60 vl>lnt(rnd(nili50+S0)«100 
70 P-lnt<rnd<l)«150+150) 
80 Pl«lnt<rnd<l)Hil50+150) 
90 c"c+l 

100 on c goto 110.150.190,230.270 
110 v«0 

120 Print: Print "Solvt for v whtn" 
130 «l-Pl««l.'P 
140 Soto 200 
150 vl«0 

160 Print "What Is vl Slwen" 
170 <Il«P*v/'Pl 
180 Soto 280 
190 P«0 

200 Print: Print "C»lcul«te v*lu» of P" 
210 <(l-Pl«vl/v 
220 9oto 280 
230 Pl-0 

240 Print: Print "Solv* for PI 91ven" 
250 <tl'P*v/vl 
260 9oto 280 
270 stop 

280 Print "v= ";v;tiib<il5; "P"";p 
290 Print "vl- ";wl;t»b(ll); "Pi=";pi 
300 InPut "Vour «nsuer";4 
310 Print "Correct v»lue";<ll 
320 Print 
330 9otoS0 
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Boy It's iMi drill 
Volunt In nilliUttr* 
Prt»»urt in can n*r 

Solvt for w when 

vl- 9980 I>1> 273 

Vour «n»wtr? 114418 
Correct v«lu* 8S21. 42837 

Uh«t U vl «lven 
V- 6480 P- 276 

vl- »1- 198 

Vour untutr? 523821 
Correct v«Iut 9296.84211 



CalcuUte vulue of P 
V- 7780 P' 8 

vl- 7400 Plm 235 

Vour »n»wer? 18494 
Correct v«lue 245.864935 

V- 7800 P- 165 

vl- 9108 PI- 8 

Vour «n»wer7 75880 
Correct vtUe 126.923877 
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Photoelectric Emissions 



When light of a short wavelength falls on a metal surface, electrons are 
emitted from the metal. According to the description of this phenomenon by 
Einstein, there is a maximum wavelength for every metal above which no 
electrons are emitted. This is called the critical wavelength of the metal. 

This program simulates a laboratory experiment in which a metal is placed 
in a vacuum and bombarded with soft X-rays. The number of electrons emit- 
ted is collected and measured with a microammeter. The program simulates 
three trials at each of nine wave lenghts. 

After each set of experimental data, the program asks if you would like an- 
other run at a higher light intensity. The reason for doing this is that some- 
times at low light intensities, not enough electrons have been emitted for 
meaningful measurements. 

You can increase the precision of the experiment by increasing the number 
of wavelengths at which it is run. This value can be changed in Statement 60; 
note that the variable L is divided into 1000 in order to express the wave- 
length in Angstroms. One Angstrom (A) equals 10"* centimeters or 10"* 
microns. 

Here are the coefficients for several metals: 



Silver .308 

Bismuth .338 

Cadmium .318 

Lead .340 

Platinum .385 



It is a rare physics laboratory in a high school or college today that has the 
experimental apparatus to run this experiment, yet with a small computer the 
equipment can be simulated. There are many other things with which you 
would not normally be able to experiment that can be simulated with a com 
puter. Things such as a nuclear power plant, a malaria epidemic, an urban 
mass transit system, and a bicycle factory. 

Can you write a simulation for a real world system? You will find sections 
on simulations in the books Computers in Mathematics and Computers in Sci- 
ence and Social Studies. These, along with articles in Creative Computing 
magazine, might be of some help in writing a simulation of your own. 



112 



5 Print chr*<14>:rtn switch loutr cue 

10 Print "Booibtrdntnt of mtttl with soft X-r«»s." 

20 InPut "Cotff lcltnt";v 

se k-lnt<l+2«rnel<l)) 

33 Print 

4fl Print "Output in mlcroamPs" 
43 Print 

30 Print " W«v» Trial Trl«V' 

35 Print "LtnSth 1 2" 

36 Print " " 

60 for 1-.42 to .23 sttP -.02, 
70 n«int( 1000^1) 

80 Print n; 

90 for J-1 to 2 

100 if I>v thtn 130 

110 i-s<lr<int<25*rnd<l)» 

120 Soto 140 

130 l>s4r(l<«k«100't'int(3S«rnci<l)>> 

140 n«int(10«i'f.S>/'10 

150 Print t«b<J#7);n; 

160 ntxt i 

170 print 

180 ntxt I 

190 Print: inPut "Incrtis* inttnsitv of light (»jn)";«* 

210 if **"''n" or »«-"N" thtn 270 

220 input "Bk wh«t factor (1 to 10)";f 

240 k-k«f 

250 Print 

260 Soto 50 

270 stop 

Bonbardntnt of nttal with soft X-r*»s. 
Coefficient? .34 
Output in nicroanPs 
Wave Trial Trial 



Ltn9th 


1 


2 


2380 


10.3 


10.2 


2500 


10.3 


10.9 


2631 


11.2 


11 


2777 


11 


11.1 


2941 


3.3 


2.8 


3125 


3.6 


4.2 


3333 


1 


4.9 


3571 


2.4 


1.7 


3846 


3.5 


1.7 



Increase intensity of U9ht (»)n>? n 
By what factor <1 to 10>? 5 

Wave Trial Trial 



LenSth 


1 


2 


2380 


50.1 


50 


2508 


50.2 


50.3 


2631 


50.1 


50.1 


2777 


50.3 


50.3 


2941 


4,7 


2.4 


3125 


4.1 


3.5 


3333 


3.9 


4.7 


3371 


2 


4.6 


3846 


i.7 


1.7 



Increase inteYisits of light (»,n)? n 
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Mutation of Moths 



This program is a simulation of the growth of a colony of pepper moths. 
The program allows a genetic mutation to be introduced in some year be- 
tween 1 and 30. The mutation can favor either dark or light colored moths. 

The program as it appears starts with a total colony of light moths; you 
could add an initial group of dark colored moths in Statement 45 as PI. 

The sample run shows a mutation which favors dark moths occurring in 
Year 3. The mutation affects about 2% of the light moths each year and 
causes them to become dark. The program displays the number of moths of 
each color over a 30-year period. 

Obviously, because of the long time period involved it would be difficult to 
carry out this experiment in school, so the computer again is of real benefit. 



S Print chi"#<14>:rtn twitch loutr cm* 
IB Print "Mut«tion of « colons of Kioth», " 
20 inPut "Ritt of mutttlon <i to 10)";m 

50 Z-P0 

£0 inPut "Envlronn«nt»l chinM occur* in vW'ix 

ee \-i 

90 cl-2 

100 inPut "Ch*n9* f«vor» U9ht or d«rk <L/ !!>";#» 
120 Print: Print "Ve«r D*rk Liflht" 

125 Print " " 

130 for t-1 to 30 
140 if t>x then 170 
150 Pl-0 
160 90to 230 

170 if t*0"D" tnd etC'd" thtn 150 

180 Pl«int(Pl+.01i|iin«p04.5) 

190 P0«int<i-Pl+.5) 

200 if PKz thtn 230 

210 Pl-z 

220 P0-0 

230 Print t;t«b<5);Pl;tab<12>;p0 
240 ntxt t 



The so-called "killer bees" that were introduced into South America in the 
mid 70's were supposed to help honey production because they were much 
more energetic than the lazy honey bees in Brazil. The idea was that they 
would mate with the existing honey bees and produce a more productive 
strain. However, they went on a rampage kilUng people and terrorizing the 
country. Then they started to migrate north. U.S. agricultural officials be- 
came alarmed that they would invade this country and bring death and 
destruction. 
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Mutition of It colora of moths. 
Rite of mutation <1 to 10)? 2 
Ettvlronneintal ch«t>9» occur* In year? 3 
Chanee favors lieht or dark <L,D)? d 



T»ar 


Dark 


Ll9ht 


4 

I 





10000 


2 





10000 


3 


200 


9800 


4 


396 


9604 


3 


588 


9412 


6 


776 


9224 


7 


960 


9040 


8 


1141 


8859 


9 


1318 


8682 


10 


1492 


8508 


11 


1662 


8338 


12 


1829 


8171 


13 


1992 


8008 


14 


2152 


7848 


15 


2309 


7691 


16 


2463 


7537 


17 


2614 


7386 


18 


2762 


7238 


* o 
19 


2907 


7093 


CO 






21 


3188 


6812 


22 


3324 


6676 


23 


3458 


6542 


24 


3589 


6411 


25 


3717 


6283 


26 


3843 


6157 


27 


3966 


6034 


28 


4087 


5913 


29 


4205 


5795 


30 


4321 


5679 



Over the years, the killer bees have started to mate with the more docile 
South American honey bees, but only at the rate of 3% per year. Assuming 
they pose a danger to the U.S. until their numbers are reduced to 15% of 
their original quantity, how many years will it take for this to occur? The 
present migration patterns will bring them to the U.S. by 1989; will they still 
be dangerous (assume that the year of introduction was 1974)? 
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Projectile Motion 



The path followed by a projectile is called its trajectory. The trajectory is 
affected to a large extent by air resistance, which makes an exact analysis of 
the motion extremely complex. We shall, however, neglect the effects of air 
resistance and assume the motion takes place in empty space. 

In the general case of projectile motion, the body (bullet, rocket, mortar, 
etc.) is given an initial velocity at some angle 6 above (or below) the 
horizontal. 

If Vq represents the initial velocity (muzzle velocity), the horizontal and 
vertical components are: 

Vox = V°s^' Voy = Vi"^ 

Since we are neglecting air resistance, the horizontal velocity component 
remains constant throughout the motion. At any time, it is: 

= Vqj^ = COS0 = constant (1) 

The vertical part of the motion is one of constant downward acceleration 
due to gravity. It is the same as for a body projected straight upward with an 
initial velocity Vq sin 6. At a time "t" after the start, the vertical velocity is: 

Vy = Voy-gt = V„sin0-gt (2) 
where "g" is the acceleration due to gravity. 




Vq sin e 
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The horizontal distance is given by: 

x = Vpxt = (V^cos9)t (3) 

and the vertical distance by: 

y=V„yt-l/2 gt^ 

= (VQSin(9)t- l/2gt^ (4) 

The time for the projectile to return to its initial elevation is found from 
Equation (4) by setting y = 0. This gives 

2 V-sin d 

' — — 

The horizontal distance when the projectile returns to its initial elevation is 
called the horizontal range. "R." Introducing the time to reach the point in 
Equation (3), we find: 

2 V^sin e cos d 

R= — (6) 

g 

Since 2sin G cos 6 = sin 26, Equation 6 becomes: 

sin 26 

R = — (7) 

g 

The horizontal range is thus proportional to the square of the initial veloc- 
ity for a given angle of elevation. Since the maximum value of sin 20 is 1, the 
maximum horizontal range, Rmax is Vg^/g. But if sin 26 = I, then 26 = 90° 
and 6 = 45°. Hence the maximum horizontal range, in the absence of air 
resistance, is attained with angle of elevation of 45°. 

From the standpoint of gunnery, what one usually wishes to know is what 
the angle of elevation should be for a given muzzle velocity Vq in order to hit 
a target whose position is known. Assuming the target and gun are at the 
same elevation and the target is at a distance R, Equation (7) may be solved 
for a 



9 = 1/2 sin' 
= 1/2 sin 



(8) 



Provided R is less than the maximum range, this equation has two solu- 
tions for values of 6 between 0° and 90°. Either of the angles gives the same 
range. Of course, time of flight and maximum height reached are both greater 
for the high angle trajectory. 
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For example, say the maximum range of our gun is 10,000 yards and the 
target is at 5,900 yards: 



Try the computer game Gunner which appears below. Use trial and error 
to try and destroy the target (see sample run). You get five chances per target. 
How many shots did it take to destroy all five targets? Did you ever fail to de- 
stroy a target in five trials? 

From the discussion above, you should realize that 45 degrees of elevation 
provides maximum range with values over or under 45° providing less range. 

The maximum range of the gun will vary between 20,000 and 60,000 yards 
and the burst radius of a shell is 100 yards. 

You can also determine the correct firing angle from sine tables, or a sUde 
rule, a scientific calculator or a computer. You should be able to destroy ev- 
ery target with just one shot. What happens when the target is very close? 
Can you always use whole angles? 

Now write a computer program to accept the maximum range of your gun 
and the range to the target and then calculate the correct firing angle. You 
will have to solve two problems to write such a program: 

1. The Basic language does not have an ARCSIN function. However, the 
following formula may help. 



e = 1/2 sin' 



,-1 



5,900 
10,000 



= 1/2 36° 

= 18°, or 90°- 18° = 72° 



sin"' X = tan' 




2.You must convert from radians to degrees. 
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5 Print chri<14>:rt« »wltch low»r c««t 
S Print chr«<147> 

18 Print "Vou *r» cMmntndlnfl « 9un crtw. Hit within 100 y«rd« of tht tir9»t." 
50 Print 

70 r-lnt<400ee«rncKl)<f20e00> 
80 Print "M*x rM9«";r;"ii»rd«" 
90 t-0 

100 Sl>0 

110 t"lnt<r«(.l+.8«rnd(l)>> 
120 f-O 
130 «oto 310 

140 Print "Nln ■ 1 dt9rtt" 
190 9oto 320 

160 Print "n«x ■ 89 d*9rtM" 
170 «oto 320 

180 Print "Ovtr bi»";»b»(t); "K«r«l»" 
190 9oto 320 

200 Print "Short b*";«b«<»)"K«rd«" 
210 Ooto 320 

220 Print "T«r9tt dntroucdl" 

230 Print i"Round* uttd" 

240 Print 

290 tl-tU» 

260 If v4 then 4S0 

270 xn*l 

280 Print 

290 Print "Mort trmnn «ctlvlt»" 
300 9oto 110 

310 Print "T»r9»t *t";t; ")*«rd»" 
320 rtn 

330 Print: Input "Eltvttlon";b 
340 If b>89 th«n 160 
390 If b<l then 140 
360 m*l 

370 If t<6 thtn 410 
380 Print chr«<147) 

390 Print t*b<12)"Boo«!"; Print t<b<9)"Booi>il": Print t«b<14) "Boom!": Print 

400 Print "EntinK 9ot vou flr»t":9oto 480 

410 b2-. e39«b : i»r«»ln<b2> : x-t-i : •-int<x> 

420 if tb*<»Km thtn 220 

430 If t>100 thtn 200 

440 9ete 188 

490 Print chr«<147);"Tot«l round* -»;»1 

460 If *1>18 thtn 480 

470 Print "Nict «hootin9! :9oto 490 

480 Prlnt:Prlnt "Oo back to Fort Sill for nort tr*inin9."; 

490 Print :lnPut "Tr* a9«ln <v,y\)"ltM 

900 if**""ii" or z*«"V" thtn 50 

910 Printiprint "OK. Rtturn to ctioP." 

920 stop 

590 tnd 
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Vou «r« conmandlnS « gun crew. Hit within 190 y«rd* of th» t«r9et. 
n«x r«n9» 33171 K»rel» 
T»r9tt «t 39743 S(»rd» 

Elfw*tlon? 28.5 
Ovtr by 4927 vtrdi 

Eltv*tlon? 25 
Ovtr bv 1269 v«rd» 

El«w«tlon7 23 
Short by 1415 y«rdt 

EUv*tlon? 23.7 
Short by 524 yard* 

El»v»tlon? 23.9 
Short by 274 yardt 

Elevation? 24 

Boon ! 

Boom! 

Boom! 

Enemy 9ot you first 

Oo bacl< to Fort SlU for wore training. 
Try a9aln <y,n>? y 



Wax ranee 58685 yards 
TarSet at 25871 sard* 

Elevation? 15 
Over by 3543 yard* 

Elevation? 12 
Short by 1941 yard* 

Elevation? 13.8 
Over by 1385 yard* 

Elevation? 12.6 
Short by 821 yard* 

Elevation? 13 
Target destroyed! 
5 Round* used 

Try aSain (B,n)? n 
OK. Return to camP. 
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9 

Potpourri 



Here are five programs that didn't seem to fit anywhere else. The first and 
the last are games, although you may come to think of the lunar landing 
simulation as a game also. One is a nifty simulation of smog, and the other 
calculates depreciation by three different methods. 
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Number Guessing Game 



Here is a computer program that plays the popular number guessing game. 
In it the computer picks a secret number between 1 and 100. You attempt to 
guess that number in as few tries as possible. 

There are many ways to go about guessing the secret number. Let some of 
your friends play this game and see what approaches they use to find the se- 
cret number. One approach is to start with a guess of 10. If this is too low, in- 
crease each guess by 10 until the computer says that a guess is too high. 
When this point is reached, start from the previous guess and increase each 
guess by 1. This is the method used to solve some of the problems in the chap- 
ter on convergence. 

Another approach is to try to bracket the number between upper and lower 
hmits and reduce the limits by steps until the number is finally found. Two of 
the convergence programs used this approach. 

Is one of these the best way? Well, these methods are not bad, particularly 
compared to starting with 1 and simply counting to 100 until the solution is 
found. But there is a better way. It is known as binary search. 

This technique involves dividing the search domain, in this case 1 to 100, in . 
half, and then in half again, and so on until the secret number is found. Play 
the game many times using different approaches. In the long run you should 
find that the binary search approach is the most efficient. 

In Line 230, the program contains the statement that "you should not need 
more than 7 guesses." Why? If the secret number was between 1 and 128, 
what is the maximum number of guesses that would be necessary to find it? 
What if the number range were 1 to 130; then what would be the maximum 
number of guesses? 

Revise the program to choose a secret number between 1 and 10,000. Now 
the upper limit is 100 times the 1 to 100 game here which requires a maxi- 
mum of seven guesses to find the secret number; how many guesses will now 
be required? 

Can you write a program in which the roles of the computer and player are 
reversed? In other words, the computer will try to guess your secret number 
between 1 and some upper limit. After each guess, you enter L for low, H for 
high, or C for correct. Can you write this program so the computer can tell if 
you are cheating, i.e., giving it inconsistent clues? 
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3 Print chrf<14):r«w switch lowtr c«ft* 

10 Print "I'vt t Mcret wnbtr bttwten 1 i 100. Tr» to 9ue»s it."; 
20 Print " I'U Slvi: KOu clut*. " 
30 Print 
40 Print 
50 c-0 

£0 n>lnt<rnd<l>liil00>'fl 
70 c-c+1 

se Print taba0);"OuMt"; 

90 inPut 9 

100 ir n-9 then 160 

110 it 9>n then 140 

120 Print "Too low." 

130 9oto 70 

140 Print "Too hi9h. " 

130 9oto 70 

1£0 Print:prlnt "Correct ir>";c;"trie»."; 

170 if c>3 then 200 

180 Print " Vou were luck»." 

190 9oto 240 

200 if c>7 then 230 

210 Print " Oood Job." 

220 9oto 240 

230 Print " Vou ihouldn't need more th»n 7 tries." 
240 Print 
250 Print 

260 inPut "PtM tnother g«ne <s«>n>";«* 
270 if «*""!(" or tl-"y" then 30 
280 stop 



I've « secret number between 1 & 100, 
Try to 9uess it. 
I'll 9ive you clues. 



Ouess? 30 

Too low. 

Ouess? 75 

Too low. 

Ouess? 85 

Too low. 

Ouess? 90 

Too low. 

Ouess? 95 
Correct in 5 tries. Oood Job. 

PlM «nother game <y,n)? i> 
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Depreciation — Three IVIethods 



This program shows how a piece of capital equipment depreciates accord- 
ing to three commonly used methods of depreciation: straight line, sum of the 
year digits, and double declining. 

The program asks for the original cost of the item, its expected life in years 
(the period of time over which it is to be depreciated), and its expected scrap 
(or sale) value at the end of that time. A table showing the annual depreci- 
ation for each of the three methods is then displayed. 



5 Print chr«<14):r«m switch lower c«»» 
le Print "Depreciation bv 3 wtthod*." 
20 Print 

30 inPut "Orl9on«t co»t";c 
50 inPut "Life <»e»rs)";l 
70 inPut "Scr«P u*lue";» 
90 Print 
100 v«c-« 
110 

120 Print "Str»i9ht line depreciation is «";eii;"Per Kear. " :print 
130 K-<<t+n/'2)#l 
140 z-l 

150 Print " Sum of Double" 
160 Print "Vear Difliti DecUnin9 

165 print " " 

170 for x»l to I 

180 dS'vKKz/u) 

190 d2«int<d2#100+.3)/'100 

200 z«z-l 

210 d3»2«c/l 

220 d3"lnt<d3*100+.5)/'100 
230 c-c-d3 

240 Print x;tab<5);d2;tab(17);d3 
250 next x 



Depreciation by 3 methods. 

OriSonal cost? 8008 
Life (years)? 8 
ScraP value? 500 

Straiflht line depreciation is * 937.5 Per year. 





Sun of 


Double 


Vear 


Digits 


Declininfl 


1 


1666.67 


2000 


2 


1458.33 


1500 


3 


1250 


1125 


4 


1041.67 


843.75 


5 


833.33 


632.81 


6 


625 


474.61 


7 


416.67 


355.96 


8 


208.33 


266.97 
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Smog Simulation 



This program is an adaptation of the smog model originally written by 
Herbert Peckham. The model assumes that vehicular traffic is the sole pro- 
ducer of smog, a somewhat poor assumption. It also assumes that most auto- 
mobile traffic occurs during the daylight hours and that traffic volume is very 
low (actually, zero) at night. The smog generated by the cars is dissipated by 
atmospheric conditions which vary depending upon sunlight, temperature, 
and weather. All of these conditions may be specified by the user. 

The model could be improved significantly by taking into account other 
sources of smog, by varying vehicular traffic according to the hour of the day, 
and by allowing daily variation of weather factors. Nevertheless, even in this 
rudimentary form it is interesting and instructive. 

A plot of the smog level is produced in Statement 400. Under some con- 
ditions, the smog level reaches a value that cannot be plotted because it is 
greater than the width of the screen (and printer). For runs with conditions 
like this, you might want to delete the plotting routine. A more elegant solu- 
tion would be to estimate the maximum value of the smog level from the in- 
put factors and calculate an appropriate plotting multiplication factor. 



10 Prlntchr*<147>;"Si«io9 SlmuUtlon" 
20 Print 
30 W 
35 Print 

40/''lnt "C«r» Ptr ro«cl iiilU":prlnt "Lot Rngtle* 200":prlnt "B*trolfc 

45 Print "Tul»» 25" 
60 Print "S«o9 Cltn ";:inPutc 
65 Print 

70 Print "Smo9 9«ntr«tlon fuctor" :prlnt "1930 fluto .005" :prlnt"1975 Ruto 

001 

75 Print "Bu» .010" 
80 Print "Smo9 Cltsi "rlnPutU 
90 Print 

100 Print "D*s(tlmt sno9 breakdown": Print "In Ptrctnt Ptr hour" 
105 Print "ClMr .02":prlnt "Clouds .01" 

110 Print "SIII09 Clt» ";:inPutrl 
120 Print 

130 Print "Nighttime breakdown" : Print "Hot .05":prlnt "Cool t«" 

140 Print "Smo9 CltK ";:inPutr2 "'^ 
150 Print 

160 Print "Dispersion <Ji Per hourV':Prlnt "Hl9h wind & rain 1.00" 
165 Print "No wind and drs 0.01" 
170 Print "Sm9 Cits* "rlnPutrS 
180 Print 

190 inPut "Table or Plot <t or P>";a* 
200 Print "Hour Si>io9 Level" 
210 r«<rl 
220 kl-ll 

230 tl«int<i:t-6V12) 

240 if tl/2"int<tl^2) then 270 
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250 ki-e 
260 r"r2 

270 »«*+kl*cilH0-r«»-r3«» 
280 If »<"0 then s"0 
290 t"t+l 

300 t3=iflt<t/'12) 

310 t2»t-t3«12+l 

320 x»lT)t<1000ilis+.5)/'1000 

330 If t3/'2«lnt<t3/2) then 330 

340 Print t2; "PM";t«b<7); :9oto 3«0 

350 Print t2; "flM";t*b<7>; 

360 if a*«"P" or »*-"P" then 390 

370 Print x 

380 9oto 210 

390 if x>17 thtn x»17 

400 Print t»b<6+x>;"#" 

410 9oto 210 



Snoo SlnuUtion 



Girt Per ro«d ml It 
Lot flnOtltt 200 

Octroit 100 
Tultt 25 

Smo9 City ? 100 

SmoO generation factor 
1990 fluto .009 
1^75 Ruto .001 
But .010 
Snog City ? .001 

Daytime tnos breakdown 
in Percent Per hour 
Clear .02 

Cloudy .01 
Smog City? .01 

Nighttime breakdown 
Hot .05 
Cool . 1 

Smog City? .1 

BitPertion <Ji Per hour) 
High wind 1 rain 1.00 
No wind and dry 0.01 
Smog City ? .05 
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T«bU or Plot <t or P>? t 



Hour 



Smo9 Ltvvl 



9 RM 


1 


le RM 


1.94 


11 RM 


2.824 


12 RM 


3.654 


1 PM 


4.433 


2 PM 


5. 169 


3 PM 


5.839 


4 PM 


6.387 


S PM 


7.117 


6 PM 


7.69 


7 PM 


8.228 


8 PM 


6.994 


9 PM 


3.943 


le PM 


5.833 


11 PM 


4.293 


12 PM 


3.631 


1 RM 


3.183 


2 RM 


2.638 


3 RM 


2.242 


4 RM 
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5 RM 


1.62 


S RM 
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7 RM 


1.17 


8 RM 


2.1 


9 AM 
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11 RM 
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12 RM 


5.294 


1 PM 


5.976 
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Lunar Lander Simulation 



This program is one of the most popular computer simulations around. It 
is available in many versions; this one is adapted from the original program 
written in 1969. 

The program represents an exact simulation of an Apollo lunar landing 
module during the final descent. This portion of the descent would normally 
be controlled by the on-board computer backed up by another computer in 
the lunar orbiter, and still another computer on Earth. However, to exercise 
your knowledge of physics and to make an interesting game, we will assume 
that all three computers have had a simultaneous malfunction. Hence, it is up 
to you to land the spacecraft safely. 

To make a soft landing, you may change the burn rate of the retro rockets 
every ten seconds. You have a choice of not firing at all (burn rate of 0) or of 
firing at a fuel rate of between 8 and 200 pounds per second. Engine ignition 
occurs at 8 pounds, hence values between 1 and 7 pounds are not possible. 
You have 16,500 pounds of fuel. This is 500 pounds more than an actual 
LEM has, which will give you a httle margin for error. When you get pro- 
ficient, change Statement 130 to N = 16000 to simulate the real thing more 
closely. The capsule weight is 33,000 pounds. 

Not that this is the way to come in, but if you did not fire the rockets at all, 
the estimated time for a free fall descent is 120 seconds to impact (and a huge 
splat). 

Good luck! 



S Print chr«(14) :rtn twltcn louer c«s* 
18 Print chr«(147);" Lun«r Undtr" 

13 Print " b» Dtvld fihl 

20 Print:Prlnt "ConPuttr iwlfunctlon - t«kt iii«m«l control of Lun«r c*P«ule. " 
30 Print 

40 Print "C«Piult ufi9ht — 32500 lb«" 
45 Print "Futl uti9ht -- 16S00 Ibt" 
50 Print 

£0 print "Stt retro burn rttt tveru 10 stcs to or a v«lut bctuetn"; 
70 Print "8 & 200 lb« Ptr ttc. " 
80 Print 

90 Print "MHt« MPH Fuel Tine Burn" 

95 Print " " 

100 «-120 
110 v-1 
120 Ki-33000 
130 n-l£500 
140 9-. 001 
150 r-1.8 

160 P-int(3£00«vf.5> 

170 If ibtCP))-! then 190 

180 P-int((3S00«V'f.5)!lil0e0)/'1000 

190 if «<10 then 210 

200 d-intU-f.S^^toto 240 

210 if «<1 then 230 
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240 Print d;Ub<8>;p;t«b(17);m-n; t*b<26>; t;Ub<34>; : inPut k 
250 If k>2e0 thtn 290 
270 If k>7 thtn 300 
280 if k-0 thtn 300 

290 Print "Ctn't do— «9«ln Ple«»e":9oto 240 
300 t*i0 

320 if m-n<.00i thtn 430 
330 if t<.001 thtn 160 
340 »-t 

330 if tOn*t#k thtn 378 
3fi0 ••<iii-n)/'k 
370 9osub 770 
380 if i<-0 thtn 630 
390 If v<«0 thtn 410 
400 if J<0 tl^n 690 
410 9osub 570 
420 9ot0 320 

430 Print "Vou r»n out of futl «t"; I; "stcond*. " 
440 »»<-v+*<lr<w*v+2«««9)>/'9 
450 v*v+9#» 
460 l-Ut 
470 u-3600iv 

480 Print "On thi »oon jit" ; lnt<n; "»tc. Ii»P«ct velocity was"; 

485 Print iTit<100#w+.5)/100; "MPh" 

490 if u>1.2 thtn 510 

500 Print "Ptrftct Undln9!":tnd 

510 if w>10 thtn 530 

320 Print "OoooonPh! That's rouSh on the shocks. Vou nttd mort Practlctl":tnd 
530 if w>60 thtn 550 

540 Print "Severt craft dana9e! Vou 're stranded until help arrives. end 

550 Print "Sorry/ there were no survivors. Vou bleu it." 

560 Print "In fact< you blasted a new lunar crater"; int<u«. 278); "feet dtepi":end 

570 iBl+s 

580 t«t-s 

590 t'i»iii-»#k 

£00 a>i 

610 v«j 

620 return 

630 if s<.005 then 470 

640 d«v+sqr<vi(iv+2«a*<9-z«k/'iii)) 

650 6»2*i./6 

660 9osub 770 

670 9osub 570 

680 9oto 630 

690 w=<l-n*e/'<z#k)V2 

780 s«i«i*v/<z*k*<up+sqr<w«w+v/z)))+.05 

710 9osub 770 

720 If l<-0 thtn 630 

730 sosub 570 

740 if J>0 then 320 

750 If v>0 then 690 

760 9oto 320 

770 <f»s«k^M 

780 J«v+9*s-z«<t#< i+q«< . 5+q«<i/'3+q*< . as+q/S) ) ) ) 

790 i aa-SfslHsr'a- v*s+z#sM *<.5+'5«<l/6+q#<l/'12*1^20>)) 
800 return 
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Lunir Under 
bv David RhI 

Computer malfunction - take nanual control of Lunar caPsute, 

Capsule wei9ht — 32308 lb> 
Fuel ueiSht — ISSee lbs 

Set retro burn rate ever» 10 sees to or a value between 8 4 208 lbs Per sec. 



Miles 


MPH 


Fuel 


Time 




Bum 


120 


3600 


1£S00 










110 


3S3£ 


1£300 


10 


? 





100 


3£72 


1£S00 


20 


■? 





90 


3708 


1£500 


30 


? 


30 


79 


3£43 


1£000 


48 


? 


30 


£9 


3381 


13500 


50 




30 


39 


3515 


13000 


£8 


? 


100 


30 


3341 


14000 


70 


? 


100 


41 


31fil 


13000 


88 


? 


100 


32 


2974 


12000 


98 


? 


100 


24 


2779 


11000 


108 


? 


200 


17 


2325 


9000 


118 


? 


200 


11 


1832 


7000 


128 


? 


200 


7.1 


1292 


5000 


138 


? 


200 


4.4 


£95 


3000 


148 


? 


200 


3.3 


30 


1000 


158 


? 


150 


Vou ran 


out of fuel at 


15£.££6££7 


seconds. 



On the noon at 372 sec. ImPact velocits* was 448.42 MPh 

Sorrs, there were no survivors. Vou blew it. 

In fact/ »ou blasted a new lunar crater 124 feet deeP ! 




Hammurabi 



Hammurabi is one of the all-time favorite computer games. On the one 
hand, it may be considered a game, but on the other it is an intriguing simula- 
tion of barter and management. 

Hammurabi is your servant as you try to manage the ancient city-state of 
Sumeria. The economy of the city-state revolves around just one thing — the 
annual crop of grain (probably soybeans). 

Each year, you must determine how many bushels of grain you wish to 
feed to your people (you'll quickly discover how much a person needs to sur- 
vive), how much you wish to use as seed in planting crops for the coming 
year, how much you wish to use for the purchase of additional land from 
your neighboring city-state, and how much you wish to put in storage. 

Of course, if you have a bad harvest or if rats overrun your grain storage 
bins, you may have to sell land in order to get enough grain to keep your peo- 
ple from starving, or to plant the land for the coming year. Unfortunately, 
disasters always seem to strike, forcing you to sell land, when the price is at 
an all-time low; but that's not any different from the real world. 

Most people start to play this game with noble ambitions. However, before 
long, they start longing for a plague to trim their growing population. Or they 
deliberately starve some people to keep things in balance (gosh, maybe these 
zero population growth people have something, after all!). 

Over the years, this game more than any other, has spawned a host of look- 
alikes, extensions, and modifications. Indeed, several manufacturers have 
taken my original with no changes whatsoever, put it in a fancy box, and 
charged a handsome price for it. Accept no imitations! Here is the original 
game (with the dialog shortened slightly) for you to run on your computer. 

If you want to experiment with changes, here are some suggestions. In the 
existing game, plagues randomly occur 15% of the time; lower this to 10% or 
5%. People now require a fixed amount of food; vary this amount slightly 
from year to year. Permit the construction of a rat-proof grain bin, but this 
must cost a fair amount. Introduce a mining industry as well as agriculture. 
How about fishing or tourism? Let your imagination run wild. Experiment! 
Have fun! 
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5 Print chr*<14>:rtm suiltch lower case 

10 Print chr*<147>;"Trv to sovern indent 8unerl« for 10 »e«rs. " 
30 dl-0:pl-e 

40 z-e:P-95:*-2800:h-3000:e-h-s 
30 i(«3:«"h/'»:l«5:^«l 
£0 d-0 

70 Print :z"r+l 

80 Print " H«nmur*bl: Ve*r";2 
90 Print d;»8t«rved/';l;"Born" 

100 P"P+i 

110 ir 4>0 then 140 
120 P»lnt<P/2> 

130 Print " Horrible Pl«<lue!" 
140 Print " PoPol»tlon -'SP 

130 Print " Vou oun";a;"icre« which Produced";sr; " bushels Per acre." 

180 Print " Rats ate" ;e; "bushels. " 

190 Print s; "bushels stored. Print 

200 If z-U then 940 

210 c"lnt<10«rnd<l)Vy-c+17 

220 Print "Land «";»;"bushels/'acre" 

230 InPut "Acres to bu»";fl 

250 If 4<0 then 810 

260 If u«<l<-s then 290 

270 «osub 730 

280 goto 230 

290 If 4-0 then 320 

300 a-a+4:s"S-K*4:c"0 

310 9oto 380 

320 Input "ficres to sell";-? 

330 If <1<0 then 810 

340 If <l<a then 370 

350 Sosub 760 

360 9oto 320 

370 a"a-4:s"s+Di(t4:c»0 

380 Print 

390 Input "Bushels for feed";q 

400 If 4<0 then 810 

410 If <l<»s then 440 

420 «osub 730 

430 9oto 390 

440 S"S-fl:c«l 

450 Print 

460 Input "Acres to setd";d:lf d"0 then 570 
470 If dC0 then 810 
480 If d<-athen 500 
490 Soto 460 

500 If lnt<d/2K-s then 530 

510 9osub 730 

520 9oto 460 

530 If d<10#P then 560 

540 Print "Not enouSh People" 

330 «oto 450 

560 s«s-lnt<d/^> 

370 Sesub 790 

580 K«c:h"d»v:e"0 

390 9osub 790 

600 If lnt<c/'2><>c/2 then 620 
620 s«s-e+h 
630 9osub 790 

640 l-lnt<c#<20ilia+s)/'P/10e+15 

650 c»lnt<q/20> 

660 <l«lnt<10*<2«rnd<l>-,3>) 
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670 lfP<c then 60 

680 d-P-cMf cl>.45«P thtn 710 

690 Pl-«z-l)«Pl+cl*100^P>/z 

700 P-c'dl-dl+d^Seto 70 

710 PrlTit:Prlnt d;"PeoPte starvtd!" 

720 Print "Vou «re Impeached. ":9oto 1030 

730 Print "Think «9»ln. Vou have Ju«t";»; "buihels" 

750 return 

760 Print "Think aSaln. Vou own onl!(";a; "acres" 
780 return 

790 c»lnt<rnd<l>«5>+l 
800 return 

810 Print "Hammurabi: I can't do that. Oet sourself another steward! " :floto 1030 

840 Print "In 10 »ears";PlJ "!i" 

850 Print "starved Per s»ear. Fl total of";dl;"dled!" 

870 l-az-p 

880 Print "Vou started with 10 acres/'Person and ended with"; 1; "acres/Person. " 
910 Print 

920 if Pl>33 then 720 
930 if K7 then 720 
940 If PD10 then 990 
950 if K9 then 990 
960 if Pl>3 then 1610 
970 if Kiethen 1810 

980 Print "Fantastic Performance! " :9oto 1030 

990 Print "Vera heavK handed! PeoPle are rebeUinS! " ^goto 1038 

1010 Print "Not bad^ but could be somewhat beter." 

1030 Print-Print "So Ions for now." 

1050 end 



Tr» to sovern ancient Suntrla for 10 Kears. 

Hammurabi: Vear 1 
Starv/tdj 5 Born 
PoPolatlon " 100 

Vou own 1000 acres which Produced 3 bushels Per acre. 

Rats ate 200 bushels. 

2880 bushels stored. 
Land * 20 bushels/acre 
Acres to buy 7 
Acres to sell? 

Bushels for feed? 1500 

Acres to seed? S00 

Hammurabi: Vear 2 
25 Starvedj 7 Bom 
PoPolatlon ■ 82 

Vou own 1000 lerts which Produced 3 bushels Per acre. 

Rats ate bushels. 

2S30 bushels stored. 
Land ■ 23 bushels/acre 
Acres to buy? 
Acres to sell? 

Bushels for feed? 1500 
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flcrt* to *ttd? see 

HAnrnurtbl^ Vtar 3 

7 St«rveel< 8 Born 
PoPoUtion ■ 93 

Vou oun leee tcrtt which Produced 2 buthtlt Ptr «ert. 

R«t» «t* buthtls. 

iSee bushtlt ttered. 
Und ■ 24 buthtlft/'«cre 
Acrts to bu»? e 
Acre* to stU? 2ee 

Bu*h*U for fMd? ISBS 

flcrM to Mcd? see 

H««nur»bi: V«»r 4 

8 St«rvtd/ fi Born 
Horrlbl* PtJi<lut! 
PoPoUtlon ■ 40 

Vou oun see ncrn which Produced 4 buthtls Ptr «crt. 

Rtt» ate e bushtlt. 

68Se buthtU ttortd. 
Und •■ 20 bu*htl*/«crt 
Rcrts to buy? 
flcrt» to mU? 

Bush* I « for fted? 2000 

ficrts to setd? 800 
Not tnouSh P top It 
flcrts to sttd? 500 
Not tnou9h PtoPlt 
ftcrts to sttd? 300 

H»iimiur»bl: Vt»r S 
St«rvtd/ 27 Born 
PoPoUtlon ■ 67 

Vou oun 800 tcrts which Productd 1 bushels Ptr «crt. 

Rtts «tt bushels. 

S000 bushels stored. 
L«nd ■ 20 bushels/tcre 
ftcres to bus*? 
Reres to sell? 20 

Bushels for feed? 2000 



Rcres to seed? S00 

Hammurabi : Vear 6 
Starved/ 9 Bom 
PoPolatlon ■ 76 

Vou own 780 acre* which Produced 2 bushels Per acre. 

Rats ate bushels. 

41S0 bushels stored. 
Land • 20 bushelSi^acre 
Rcres to buy? 
Rcres to sell? 80 

Bushels for feed? 2300 
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flcr»» to seed? 508 

Hiiiifiiurjibl : Ve»r 7 
Starved < 8 Born 
PoPolntlon ■ 84 

Vou own 700 «cres which Produced 2 bushels Per »cre. 

R«t* it* bushels. 

4000 bushels stored. 
L«nd • 22 bushel s/'acre 
Acres to buM? 
ficres to sell? 

Bushels for feed? 2000 

Rcres to seed? 880 
Acres to seed? 30 

H«i»«ur»bi: Ve«r 8 
Starved < 6 Born 
PoPoUtion " 90 

Vou own 700 acres which Produced 3 bushels Per acre. 

Rats ate bushels. 

2075 bushels stored. 
Land ■ 21 bushel si'acre 
Acres to buv? 
Acres to sell? 300 

Bushels for feed? 3000 

Rcres to seed? 300 

Hammurabi : Vear 9 
Starved) 5 Born 
PoPolatlon « 95 

Vou own 400 acres which Produced 1 bushels Per acre. 

Rats ate bushels. 

5525 bushels stored. 
Land » 18 bushel s/'acre 
Acres to bu»? 
Acres to sell? 400 
Think a9aln. Vou own on Is 400 acres 
ficres to sell? 40 

Bushels for feed? 2000 

Rcres to seed? 500 
Acres to seed? 300 

Hammurabi : Vear 10 
Starved) 7 Born 
PoPolatlon " 102 

Vou own 360 acres which Produced 2 bushels Per acre. 

Rats ate bushels. 

4695 bushels stored. 
Land ■ 19 bushels/acre 
Acres to bus? 100 

Bushels for feed? 

Acres to seed? 
102 People starved! 
Vou are Impeached. 

So lon9 for now. 
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The magazine referred to in the text is: 

• Creative Computing. This is the leading magazine of software and appUca- 
tions for all small computers. It carries articles, tutorials, how-to applications, 
and extensive in-depth evaluations. 

Books referred to in the text include: 

• Computers in Mathematics: A Sourcebook of Ideas. Hundreds of classroom- 
tested ideas for using computers to learn about mathematics. 

• Computers in Science and Social Studies. Scores of simulation programs in 
biology, ecology, physics and management of real world systems. 

All of these books and magazines are available from Creative Computing. 
Write or call for the current price: 

Creative Computing 

39 E. Hanover Ave. 

Morris Plains, NJ 07950 

(800) 631-8112 

In NJ (201) 540-0445 
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